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Volume  I— DESCRIPTION  AND  INSTRUCTIONS  FOR  USE  OK  THE  AUTOMATED  PROCEDURE 
Volume  n— LISTING  AND  GLOSSARY  FOR  PROGRAM  SLOPimP 


Computer  progrfuns 
Mapping 

SLOPhUAP  (Computer  program) 


Slopes 

Topographic  map.a 


lACT  fCaml^um  as  mw9tm  aM*  tt  Af  AiocA  mm^m) 

An  automated  procedure  for  constructiiig  slope  maps  was  developed  and  dem^ 
ited  in  this  study.  The  procedure  consists  of  tliree  sequential  parts; 
iput  of  surface  elevation  values  by  means  of  a matrix  of  elevation  values 
■ed  to  as  an  elevation  grid  array,  (2)  calculation  of  slope  magnitudes 
rections  by  the  computer  program  SLOPkWAl',  and  (3)  construction  ot  slope 
ising  various  SLOPEIMAP  output  products. 
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20.  ABi'TKACT  (Continued). 

The  form  of  Input  elevation  data  selected  for  use  by  SliOPOlAI'  to  be  pre- 
pared In  advance  of  program  execution  Is  an  orthogonal  two-dimensional  matrix, 
each  element  of  which  Is  the  elevation  of  the  tojiogrnphlc  surface  at  that 
r.ntrlx  (or  grid)  point.  At  run  time,  SLCiriMAl’  Inputs  these  elevation  values 
nr.  i calculates  the  slope  magnitude  and  direction  of  each  input  grid  point  by 
employing  approximation  methods  that  use  the  elevation  values  of  the  grid 
point  in  question  and  Its  nearest  and  next  nearest  neighbors.  Then,  depending 
on  the  selected  user  options,  the  calculations  are  recorded  or  displayed  In 
ne  or  more  of  the  following  forms:  (l)  printer  swath  dumps,  (2)  magnetic 
tapes,  (3)  punched  cards,  and/or  [^)  drum  or  cathode  ray  tube  (CRT)  graphic 
; lots  that  delineate  slope  classes  selected  by  the  user. 

Volume  I of  this  two-volume  report  contains  detailed  Instructions  de- 
scribing the  execution  of  SLOPOIAP  and  two  examples  demonstrating  the  slope 
nap  construction  procedure.  Volume  II  provides  a computer  listing  of  the 
program  SLOPEMAP  and  a glossary  of  program  variables.  SUOPDIAI’  wan  coded 
with  Fortran  IV  to  ensure  future  implementation  ease  on  other  computer 
systems. 
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PREFACE 


The  study  reported  herein  was  conducted  under  two  projects:  (a) 
Department  of  the  Arriy  (DA)  Project  **A152121A896,  "Environmental  Quality 
for  Construction  and  Operation  of  Military  Facilities,"  Task  01,  "En- 
vironmental Quality  Management  for  Military  Facilities,"  Work  Unit  006, 
"Methodology  for  Characterization  of  Military  Installations  Environ- 
mental Baseline,"  sponsored  by  the  Directorate  of  Military  Construction, 
Office,  Chief  of  Engineers  (OCE);  and  (b)  DA  Project  IE865803MT30, 
"Research  and  Development  Computer  and  Information  Science  Support," 
sponsored  by  the  Engineer  Information  and  Data  Systems  Office,  OCE, 
under  the  "Improved  Data  Effectiveness  and  Availability"  program. 

.*>lr.  Vincent  Gottschalk  was  the  OCE  Technical  Monitor. 

The  study  was  conducted  during  the  period  1 June  1975  to  30  June 
1976  by  personnel  of  the  Environmental  Simulation  Branch  (ESB),  Environ- 
mental Systems  Division  (ESD),  Mobility  and  Environmental  Systems  Lab- 
oratory (MEGL),  U.  S.  Army  Engineer  Waterways  Experiment  Station  (WES). 
The  work  was  under  the  direct  supervision  of  Dr.  V.  E.  LaCarde,  ESD, 
and  Mr.  J.  K.  Stoll,  Chief,  ESB,  and  under  the  general  supervision  of 
Messrs.  B.  0.  Benn,  Chief,  ESD,  and  W.  G.  Shockley,  Chief,  MESL. 

Drs.  H.  Struve,  ESB,  and  LaGarde  were  responsible  for  the  slope  map 
methodology  and  associated  computer  programs,  and  Dr.  Struve  prepared 
the  report. 

COL  G.  H.  Hilt,  CE,  and  COL  J.  L.  Cannon,  CE,  were  Directors  of 
the  WES  during  the  study  and  report  preparation.  Mr.  F.  R.  Brown  was 
Technical  Director. 
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COHTEIITS 

VOLUME  II 

LISTING  AND  GLOSSARY  FX)R  PROGRAM  SLOPQ4AP 
(published  under  separate  cover) 


CONVERSION  FACTORS,  U.  S.  CUSTOMARY  TO  METRIC  (Sl)  AND  METRIC  (SI) 
TO  U.  S.  CUSTOMARY  UNITS  OF  MEASURSMIirr 


Units  of  measurement  used  in  this  report  can  be  converted  as  follows: 


Multicl 


To  Obtain 


AN  AtITOMATKD  PROCEDURE  FOR  SLOPK  MAP  COWSTRUCTIOH 


VOLUME  I:  DESCRIPTION  AND  INSTRUCTIONS  FOR  USE  OF 
THE  AUTOMATED  PROCEDURE 

PART  I:  INTRODUCTION 

BackitromtJ 

1.  It  is  the  policy  of  the  Department  of  the  Army  to  implement 
rigorously  the  National  Environmental  Policy  Act  of  1969  and  its  deriva- 
tive statutes  as  outlined  by  Executive  Order  and  other  Federal  and  state 
environmental  statutes.  This  policy  requires  that  the  Army  must  make  a 
deliberate  and  conscientious  effort  to  carry  out  its  mission  so  as  to: 

Preserve  unique  and  important  ecological,  aesthetic,  and 
cultural  components  of  our  national  heritage. 

Conserve  and  use  wisely  the  natural  resources  of  our 
nation  for  the  benefit  of  present  and  future  generations. 

£•  Restore,  maintain,  and  enhance  the  natural  and  man-made 
environment  in  terms  of  its  productivity,  variety,  spa- 
ciousness, beauty,  and  other  measures  of  quality. 

Create  new  opportunities  for  the  American  people  to  use 
and  enjoy  their  environment. 

2.  Reflections  on  the  many  and  diverse  activities  carried  out  by 
the  Array  that  could  adversely  impact  on  items  a through  d above  make  it 
obvious  that  Army  planners  (at  all  levels)  need  readily  accessible  data 
on  a large  variety  of  environmental  elements  including  topography, 
soils,  geology,  vegetation,  climatology,  surface  water,  groundwater, 
land  use,  and  others.  Further,  for  the  most  part,  these  baseline  data 
must  be  in  a form  compatible  with  modern  computer  facilities  because  of 
the  vastly  expanded  requirement  for  data  analysis  and  display.  The 
managers  of  the  Army's  research  project  "Environmental  Quality  for 
Construction  and  Operation  of  Military  Facilities"  have  adopted  a con- 
cept for  environmental  data  analyses  that  includes  the  extensive  use 

of  computerized  mathematical  prediction  models;  therefore,  a prime  cri- 
terion for  the  selection  of  an  element  to  be  incorporated  in  a baseline 


6 


duta  base  Is  its  rt* qui rerntfrit  for  th«‘  |)r«*dIctIon  of  envlronwiital  effects 
caused  by  an  Army  activity. 

3.  This  report  describes  a new  c'miputerized  procedure  for  napping 
slope,  a topoKraphic  baseline  element.  Uata  on  terrain  slop*?  is  needed 
for  a wide  variety  of  uses  including  procedures  for  evaluatln#^  the  Im- 
pact on  the  envirorunent  of  military  activities,  such  as  construct J on , 
field  maneuvers,  and  maintenorice  (preparation  ol'  environmental  impact 
assessments  (EIA'o)  and  statements  (EIS's)).  Commonly,  plruis  for  the 
mitij^ation  of  adverse  effects  must  t^ike  Into  account  the  effect  of 
terrain  slope  on  such  things  as  soil  erosion  (see  Reference  l),  surface 
runoff,  land  revegetation,  wildlife  habitat  enhancement  and  preserva- 
tion, equipment  operation  for  Imid  restoration,  and  facilities  siting. 
Further,  data  on  terrain  slopes  are  frequently  needed  in  general  terrain 
intelligence  products,  in  correcting  satellite  and  aircraft  scanner  data, 
and  in  designing  and  evaluating  most  ground  weapons  systems  and  all 
military  ground-contact  vehicles. 

I4.  The  manual  techniques  currently  used  for  mapping  slopes  from 
topographic  maps  are  slow  and  yield  somewhat  inconsistent  results.  Even 
under  optimum  circumstances,  there  exists  a probability  of  error  in 
manually  constructed  slope  maps  simply  because  of  the  many  subjective 
decisions  required  in  the  procedure.  The  resulting  cost  of  constructing 
slope  maps  manually,  as  well  as  the  inherent  risk  of  error  in  the  final 
product,  indicates  clearly  the  need  to  develop  faster  and  more  reliable 
techniques  for  producing  slope  maps. 

5.  The  Defense  Mapping  Agency  Topographic  Center  (DMATC)  has 

2 

developed  a system  of  computer  programs  that  uses  a grid  network  of 
ground-surface  elevations  and  converts  these  elevations  into  a corre- 
sponding network  of  slope  values.  The  DMATC  procedure  is  based  on  the 
least  squares  fit  of  a plane  in  the  direction  of  maximum  inclination  and 
does  not  yield  slope  data  of  an  accuracy  that  is  required  in  many  prob- 
lems dealing  with  the  impact  of  military  activities  on  the  environment, 
llie  least-squares-fit  method  (designated  plane  algorithm  in  this  report) 
was  evaluated  and  compared  to  two  other  calculation  methods.  The  re- 
sults of  this  comparative  evaluation  are  presented  in  Part  II. 
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6.  The  Initial  steps  in  develop I n#^  an  Improved  method  of  automat- 
ing slope  calculations  and  slope  map  construction  had  already  been  com- 
pleted by  personnel  at  the  U.  3.  Army  Engineer  Waterways  Experiment  Sta- 
tion (WES)  prior  to  the  initiation  of  this  project.  These  steps  include 
the  procedures  necessary  for  representing  a topographic  surface  in  a 
form  tiiat  is  compatible  with  digital  computers.  The  form  currently  used 
at  the  WKj  is  an  orthogonal  two-dimensional  matrix,  each  element  of 
which  is  the  elevation  of  the  topographic  surface  at  that  matrix  (or 
grid)  point.  The  two  dimensions  of  the  matrix  are  designated  X and 

Y , and  the  elevation  of  the  surface  is  designated  Z . Since  the  X 
and  Y values  are  understood  from  the  position  in  the  matrix,  only  the 
Z values  are  actually  recorded.  The  matrix,  frequently  referred  to  as 
an  elevation  grid  array,  can  be  produced  from  source  data  obtained  from 
field  survey  data,  aerial  photography,  or  topographic  maps. 

7.  A convenient  source  of  elevation  data  frequently  used  at  the 
WES  is  a topographic  map  of  the  area  of  interest.  The  semiautomated 
procedure  for  generating  an  elevation  grid  array  using  topographic  maps 
requires  that  one  manually  move  a cursor  along  each  contour  line  of  the 
topographic  map  and  record  X , Y , and  Z coordinates.  This  process, 
called  digitizing,  produces  a computer  magnetic  tape  of  stored  X , Y , 
and  Z values.  The  tape  can  then  be  submitted  as  input  data  to  a WES 
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computer  program  to  generate  the  final  elevation  grid  array.  Automated 
techniques  for  digitizing  the  contour  lines  by  optical  scanning  processes 
are  now  being  developed  at  the  WES  and  should,  when  completed,  signifi- 
cantly reduce  the  time  and  cost  of  producing  elevation  grid  arrays. 

Purpose 

8.  The  purpose  of  the  study  was  to  develop  and  demonstrate  the 
use  of  an  automated  procedure  for  constructing  slope  maps  from  a digital 
topographic  data  base,  i.e.  from  a terrain  surface  defined  by  an  eleva- 
tion grid  array  (see  pareigraph  6),  so  that  personnel  at  Army  installa- 
tions could  produce  cost-effective  baseline  slope  data  (i.e.  slope  maps) 
to  support  a variety  of  military  activities  requiring  preparation  of 
Environmental  Impact  Assessments  and  Statements  (EIA* s/EIS*s ) . 
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Approach  and  Scope- 


9.  The  approach  was  to  initially  select  several  methods  of  cal- 
culating slope  values  from  an  elevation  grid  array.  Next,  it  was  nec- 
essary to  develop  a computer  algorithm  for  each  method  and  determine 
which  algorithm  would  best  describe  slopes  on  most  three-dimensional 
terrain  surfaces.  The  selection  of  the  best  algorithm  was  based  on 
comparing  exact  slope  values  and  algorithm  slope  values  computed  from 
three  synthetic  elevation  grid  arrays  derived  from  three  different 
mathematically  defined  surfaces.  The  equations  representing  each  sur- 
face were  precisely  defined  and  continuous.  This  me€Uit  that  there  w«re 
no  sudden  Jumps  in  the  equations,  and  all  the  points  on  the  curves  were 
differentiable.  Therefore,  at  any  point  on  the  surfaces,  the  derivative 
(slope)  could  be  found.  The  deriva^'ves  derived  from  these  equations 
were  the  "exact  slopes."  C^ce  constructed,  the  synthetic  grid  arrays 
were  used  as  input  arrays  to  each  slope  algorithm.  The  slopes  computed 
with  each  algorithm  were  then  compared  with  the  exact  slopes  on  a grid 
point  by  grid  point  basis,  and  the  algorithm  that  produced  slope  values 
in  agreement  most  frequently  with  the  exact  values  was  selected. 

10.  After  the  optimum  algorithm  had  been  selected,  an  automated 
procedure  was  developed  to  output  the  slope  magnitudes  in  numerical  and 
classed  forms  on  a high-speed  printer,  magnetic  tape,  and  computer  piuich 
cards.  Subroutines  also  were  devt..oped  to  output  specified  slope  class 
ranges  graphically  on  drum  or  cathode  ray  tube  (CFT)  plots  and  to  output 
slope  vectors  in  component  form  on  magnetic  tape. 

11.  The  operation  and  demonstration  of  the  automated  procedure  is 
described  in  detail  in  Volume  I of  this  two-volume  report.  Two  widely 
different  terrain  surfaces  were  selected  for  demonstrating  the  pro- 
cedure. The  computer-constructed  slope  maps  for  the  two  sites  were 
compared  with  manually  constructed  slope  maps  in  terms  of  accuracy  and 
cost.  Volume  II  contains  a FORTRAN  IV  listing  of  the  computer  program 
and  a glossary  of  program  variables. 


PART  II:  DEVELOPING,  TESTING,  AND  EVALUATING  COMPUTER 
ALGORITHMS  FOR  CALCUUTING  SLOPE 


12.  To  accomplish  the  objective  of  this  study,  it  was  first 
necessary  to  select  several  different  methods  for  calculating  terrain 
slopes;  three  methods  were  selected,  and  a computer  algorithm  was 
developed  for  each.  The  slopes  calculated  with  the  algorithms  were  then 
compared  with  exact  slope  values  determined  by  differential  calculus  for 
three  mathematically  defined  surfaces,  and  the  algorithm  that  provided 
the  optimum  trade-off  between  accuracy  and  cost  was  selected  for  lise  in 
the  automated  procedure  for  slope  map  construction. 

Theoretical  Considerations 


13.  The  theoretical  considerations  relevant  to  slope  calculation 
are  briefly  discussed  in  this  section  to  provide  a common  mathemati- 
cal framework  from  which  all  the  slope  algorithms  will  be  developed. 

lU.  Points  on  a three-dimensional  terrain  surface  can  be  described 
by  an  equation  of  the  form  z = f(x,y)  , where  z is  defined  as  the 
perpendicular  distance  from  the  terrain  surface  to  a point  P(x,y) 
lying  in  the  X-Y  reference  plane  of  a Cartesian  coordinate  system. 

The  perpendicular  distance  z (Figure  l)  is  then  the  elevation  of  the 
terrain  surface  at  the  point  P(x,y,z)  with  respect  to  the  X-Y 
reference  plane. 

15.  If  the  equation  of  the  terrain  surface  is  expressed  in  its 
more  general  form,  i.e. , F(x,y,z)  = c , where  the  constant  c is 
set  equal  to  zero  so  that  F(x,y,z)  i z - f(x,y)  , it  can  be  shown 
that  the  gradient  of  F(x,y,z)  is  a vector,  N , normal  to  the  surface 
F(x,y,z)  at  the  point  P(x,y,z)  . The  vector  can  then  be  expressed 
as 


= 


^F  = - 


3f 

3y 


J + k 


(1) 
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NOTE  THE  PEROENOICUL  AR  DISTANCE  I (ORIGINATING  AT  THE 

POINT  P («,  yl  ON  THE  X-V  REFERENCE  PLANE)  REPRESENTS 
the  elevation  of  the  point  P I»,  y,  II  ON  THE  SURFACE 
F l«,  y,  II.  THE  SURFACE  NORMAL  VECTOR  N IS  SHOWN 
PERPENDICULAR  TO  THE  SURFACE  F IM,  y,  I)  AT  THE  POINT 
P l»,  y,  II  (ORIENTATION  angles  OF  fi  ARE  SHOWN  IN 
FIGURE  21 


Figure  1.  Theoretical  representation  of  a three-dimensional 

terrain  surface 

where 

^ = the  Cartesian  vector  differential  operator  del 

i,J,k  = the  orthogonal  unit  vectors  directed  along  the  X,  Y,  and  Z 
axes  of  the  Cartesian  coordinate  system 

Figure  2 shows  an  example  of  normal  to  the  surface,  together  with 
its  orientation  angles  i^  and  0 , which  are  defined  as  the  zenith  and 
azimuth  angles  of  , respectively.  Since  the  terrain  slope  and  as- 
pect angles,  defined  by  the  magnitude  and  direction  of  tilt  from  the 
horizontal  of  a plane  perpendicular  to  , are  equivalent  to  the 
zenith  and  azimuth  angles,  respectively,  no  distinction  is  made  in  the 
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terminology  for  these  angles.  For  brevity,  N is  referred  to  in  this 
study  as  a "slope  vector."  The  slope  and  aspect  angles  can  now  be 
calculated  from  the  components  of  N , and  it  follows  (Figure  2)  that 


tan  0 


V3f/3y/ 


Equation  3 was  derived  with  the  Y axis  assumed  to  be  directed  toward 
true  north  so  that  the  conventional  implied  meaning  of  aspect  angle 
would  be  maintained. 

Development  of  Algorithms 

16.  The  development  of  each  computer  algorithm  centers  about 
methods  of  calculating  the  expressions  for  the  tangent  of  the  slope  and 
aspect  angles  in  Equations  2 and  3.  Before  this  can  be  achieved,  how- 
ever, an  analytical  expression  must  first  be  determined  for  f(x,y) 
that  approximates  the  terrain  surface,  and  second,  the  partial  deriva- 
tives of  f(x,y)  with  respect  to  x and  y must  be  obtained.  In 
this  study,  the  approximations  of  f(x,y)  were  restricted  by  the  for- 
mat of  the  elevation  input  data.  The  actual  form  of  the  input  data 
will  be  discussed  more  completely  in  Part  III,  but  the  essentials  to  be 
understood  at  this  point  are  that  the  data  were  formatted  prior  to 
execution  of  the  automated  procedure  into  an  elevation  grid  array  (see 
paragraph  6)  with  the  same  fixed  grid  spacing  in  both  the  x and  y 
directions  and  with  an  elevation  value  assigned  to  each  grid  location. 
Thus,  the  development  of  the  three  slope  algorithms  discussed  in  the 
following  paragraphs  necessarily  conforms  to  these  requirements  for 
input  of  the  elevation  data. 

Vector  algorithm 

17.  The  vector  algorithm  was  selected  for  development  first 
because  of  its  conceptual  and  programming  simplicity.  The  eLLgorithm 
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F(x,y,z)  = z - f(x,y))  and  their  trigonometric  relations  with  the  components 


established,  as  a local  origin,  the  grid  p)oint  whose  slope  and  aspect 
angles  were  to  be  evaluated  and  performed  a search  around  the  grid 
point  using  the  four  nearest  neighbors  and  the  four  next-nearest  neigh- 
bors as  vector  termination  points  (Figure  3).  The  slope  and  aspect 
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Figure  3.  Reference  plane  configviration  for 
vector  algorithm 

angles  were  then  calculated  for  each  vector  position  using  as  partial 
derivatives  in  Equations  2 and  3 the  analytical  expressions  given  in 
Table  1.  The  maximum  slope  angle  determined  from  these  sets  of  deriva- 
tives was  selected  to  represent  the  slope  of  the  given  grid  point.  If 
the  grid  point  was  a map  boundary  point,  the  maximum  slope  angle  was 
selected  from  among  vectors  only  within  and  along  the  map  boundary. 

Plane  algorithm 

l8.  The  scheme  followed  for  the  plane  algorithm  was  to  select  the 
plane  with  the  greatest  slope  value  from  among  four  surrounding  planes. 
The  planes  were  defined  by  the  adjacent  nearest-neighbor  grid  points 
(Figure  and  the  grid  point  being  evaluated.  The  slope  and  aspect 
angles  were  calculated  from  the  expressions  for  the  partial  derivatives 


^ A 

Figure  i*.  Reference  plane  configuration  for 
plane  algorithm 

of  f{x,y)  given  in  Table  1.  Again,  if  the  grid  point  to  be  evaluated 
was  a map  boundary  p>oint,  the  maximum  slope  was  selected  from  among 
planes  only  within  the  map  boundary. 

Three-dimensional  surface  algorithm 

19.  The  surface  algorithm  was  used  to  calculate  the  maximum  slope 
from  two  fitted  three-dimensional  surfaces  S and  S’  . The  nearest- 
neighbor  grid  points  and  the  grid  point  to  be  evaluated  determined  the 

S surface  (Figure  5a),  and  the  next-nearest-neighbor  grid  points  and 
the  grid  point  to  be  evaluated  determined  the  S'  surface  (Figure  5b). 
The  slope  value  assigned  to  the  central  grid  point  was  calculated  from 
the  maximum  partial  derivatives  of  f(x,y)  derived  from  both  surfaces 
(Table  l).  The  maximum  partiaJ  derivatives  did  not  necessarily  come 
from  the  same  surface.  Thus,  two  independent  opportunities  were  avail- 
able for  determining  maximum  3f/3x  and  maximum  3f/3y  by  the  surface 
algorithm. 

20.  The  fitting  procedure  for  each  surface  was  accomplished  by  a 
second-order  Taylor  series  expansion.  The  central  grid  point  was 
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z 


b.  Next-nearest-neighbor 
surface  S' 


Figure  5.  Surfaces  S and  S'  as  defined  for  the  three  dimensional- 

surface  algorithm 


defined  as  a local  origin,  and  the  expansion  was  performed  about  this 
origin.  The  resulting  surface  equations  were  developed  as  described  in 
the  following  paragraphs. 

21.  Let  the  equation  z = f(x,y)  represent  the  actual  terrain 
surface,  where  z is  the  elevation  of  the  grid  point  with  respect  to 
the  X-Y  Cartesian  reference  plane.  Expanding  f(x,y)  in  a Taylor 
series  to  the  second  order  about  the  local  origin  yields 


2 2 

f(x,y)  = Zq  + ax  + by  + cx  + dxy  + ey 


(M 


where 


Zq  = the  elevation  of  the  grid  point  at  the  local  origin 
_ 3fl 


3y 


c i 1/2 


sff 

3x^ 
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The  symbol  |q  implies  that  the  value  of  each  partial  derivative  is  to 
be  evaluated  at  the  point  x = y = 0 . 

22.  Now  let  two  coordinate  systems  be  defined  whose  origins  are 
common  at  the  grid  point  of  interest.  First,  the  unprimed  system  (S), 
(Figure  6a)  is  defined  by  the  nearest-neighbor  grid  points,  and  second, 
the  primed  system  (S')  (Figure  6b)  is  defined  by  the  next-nearest- 
neighbor  grid  points.  An  equation  similar  to  Equation  1+  can  be  written 
for  the  primed  system: 


f(x',y')  = + a'x*  + b'y'  + c'x'^  + d'x'y'  + e'y'^ 


(5) 


Choose  now  a grid  spacing  of  D , and  apply  Equation  U to  each  of  the 
four  nearest-neighbor  elevation  values.  Four  equations  result: 


2 

= Zq  + aD  + cD 

2 

z^  = Zq  + bD  + eD 

2 

z^  = Zq  - aD  + cD 

2 

z.^  = Zq  - bD  + eD 


(6) 

(7) 

(8) 
(9) 


where  the  values  a , b , c , and  e are  defined  in  paragraph  20. 

Since  the  slope  calcination  requires  Just  the  values  of  3f/3x  and 
3f/3y  , values  for  a and  b can  be  derived  by  treating  Equations  6-9 
as  simultaneous  expressions,  then 


and 


a 


2D 


(10) 
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a.  Nearest-neighbor  or  unprimed  system 


Fig\ire  6,  Taylor  series  expansion  coordinate  systems. 

tions  conform  to  a right-handed  coordinate  system) 


b 


(11) 


23.  A similar  set  of  four  equations  can  be  developed  from  Equa- 
tion 5 and  the  next-nearest-neighbor  elevation  values.  The  expansion 
coefficients  a',  b' , c',  and  e'  can  be  determined,  but  as  before, 
only  the  values  of  a'  and  b'  are  needed.  They  are  expressed  as 

(z„  - z^) 

a-=^ ^ (12) 

2'f2  D 


2*^  D 


2h.  Before  the  maximum  derivatives  with  respect  to  x and  y 
could  be  determined,  it  was  necessary  to  express  Equations  12  and  13 

in  terms  of  the  unprimed  coordinates.  This  was  accomplished  by  trans- 

* 

forming  these  equations  by  a simple  U5-deg  rotation. 

25.  A summary  of  the  partial  derivatives  expressed  in  Equa- 
tions 10  and  11  and  the  transforms  of  Equations  12  and  13  are  presented 
in  Table  1.  If  the  grid  point  to  be  evaluated  was  a map  boundary  point, 
it  was  necessary  to  generate  pseudo-elevation  values  before  the  partial 
derivatives  in  Table  1 could  be  calculated.  These  pseudo-elevation 
values  were  obtained  by  reflecting  inner  grid  elevation  values  through 
the  map  boundary  point  to  corresponding  pseudo-grid  locations  outside 
the  boundary. 

Testing  of  Algorithms 
Generation  of  terrain  surfaces 

26.  To  compare  the  three  algorithms  in  a consistent  and  unbiased 
way,  it  was  necessary  to  generate  controlled  elevation  grid  arrays  that 


converting  U.  S.  customewy  units  of  measure- 
.s  and  metric  (SI)  units  to  U.  S.  customary 
e 5. 
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could  be  easily  modified  to  represent  a variety  of  terrain  features  and 
at  the  same  time  have  the  slope  and  aspect  angles  knovm  precisely  at  any 
given  location.  To  do  this,  three  mathematical  surfaces  were  developed, 
which  had  the  form  z = f(x,y)  , where  z was  the  elevation  value  and 
f(x,y)  was  a completely  continuous  function  in  x and  y . Continuity 
implied  the  existence  of  precisely  determinable  slope  values  for  any 
given  coordinate  location  of  x and  y . The  three  surfaces  generated 
were  named  cuspfunction,  cosfunction,  and  deltafunction  after  their 
analytic  functional  forms  (Figure  J), 

27.  Cuspfunction  surface.  The  cuspfunction  surface  was  generated 
by  the  exponential  form 


z 


Be 


-Ar 


(lU) 


where  A and  B are  adjustable  parameters.  The  scalar  r is  assigned 
the  magnitude  of  the  radius  vector  r = (x  - + (y  - yQ)J  , where 

Xq  and  y^  are  convenient  grid  coordinates  for  the  origin  of  the 
radius  vector  r . Fig\ire  7a  shows  a three-dimensional  perspective 
plot  of  this  surface. 

28.  Cosfunction  surface.  The  cosfunction  surface  was  generated 
by  a superposition  of  cosine  functions  and  is  given  by 


(15) 


where 

A = amplitude  in  the  x direction 
B = amplitude  in  the  y direction 
T^  = periodicity  factor  in  units  of  x 
T^  = periodicity  factor  in  units  of  y 
Figure  7b  presents  a perspective  plot  of  this  surface. 

29.  Deltafunction  surface.  The  deltafunction  surface  was  ex- 
pressed in  its  analytical  form  as 


A sln(2ffr/T)  ^ 

r 0 


(16) 
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a.  Perspective  plot  of  the  generated  cuspfur»ction  surface 
Figure  T.  Mathematical  surfaces  generated  to  test  algorithms  (sheet  1 of  3) 


Perspective  plot  of  the  generated  deltafunction  surface 


where  the  function  is  characterized  by  its  amplitude  A , periodicity 
T , and  the  scalar  magnitude  of  the  radius  vector  r = (x  - XQ)i 


+ (y 


yQ)j  . The  reference  plane  determined  by  z^ 


was  set  so  that 
Figure  Tc  shows  a perspective  plot  of 


all  z values  would  be  positive, 
this  surface. 

30.  Equations  II4-16  are  continuous  and  differentiable;  therefore, 
exact  analytical  expressions  for  the  slope  and  aspect  angles  could  be 
obtained.  With  this  capability,  the  slope  algorithms  could  be  compared 
and  evaluated. 

Testing  the  slope  algorithms 

31.  An  elevation  grid  array  was  generated  from  the  equations  for 
the  three  test  surfaces.  For  testing  purposes,  slope  class  ranges  were 
set  at  10-deg  increments  beginning  with  0 (i.e.  0 to  <10  deg,  10  to 
<20  deg,  20  to  <30  deg,  etc.),  and  an  array  of  classed  slope  values  was 
generated  from  the  elevation  data.  Slope  classes  were  delineated  by  the 
drum  plot  subroutine  (see  paragraphs  56  and  57).  Figures  8,  9,  and  10 
show  the  boundaries  of  the  delineated  slope  classes  for  the  cuspfunction, 
cosfunction,  and  deltafunction  surfaces,  respectively,  together  with 
boundaries  of  the  slope  classes  determined  by  differential  calculus 
(exact  values).  The  classes  are  labeled  according  to  the  ranges  in  the 
following  tabulation: 


Class 

A 

B 

C 

D 


Range,  deg 
0 to  <10 
10  to  <20 
20  to  <30 
30  to  <U0 


Class 

E 

F 

G 

H 


Range,  deg 
ho  to  <50 
50  to  <60 
60  to  <70 
70  to  <80 


Evaluation  of  Algorithms 

32.  Although  the  boundaries  between  slope  classes  for  all  algo- 
ritluns  displayed  the  same  trends  and  lay  very  close  to  the  boimdaries 
determined  by  differential  calculus  (exact  values),  the  boundaries 
generated  by  the  surface  algorithm  clearly  coincided  more  closely  with 
those  determined  by  differential  calculus.  Not  only  did  the  boundaries 
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a.  Determined  by  the  vector 
algorithm 


b.  Determined  by  the  plane 
algorithm 


c.  Detcimined  by  the  surface 
algorithm 


d.  Determined  by  differential 
calculus 


Fig\ire  8.  Comparison  of  slope  class  delineation  boundaries  for 
the  surface  cuspfunction 


b.  Determined  by  the  plane 
algorithm 


a.  Determined  by  the  vector 
algorithm 


d.  Determined  by  differential 
calculus 


c.  Determined  by  the  surface 
algorithm 


Figure  9.  Comparison  of  slope  class  delineation  boundaries  for 
the  surface  cosfunction 


a.  Determined  by  the  vector 
algorithm 
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b.  Determined  by  the  plane 
algorithm 
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c.  Determined  by  the  surface 
algorithm 


d.  Determined  by  differential 
calculus 


Figure  10.  Comparison  of  slope  class  delineation  boundaries  for 
the  surface  deltafunction 
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coincide  more  closely,  but  also  the  slope  values  calculated  with  the 
surface  algoritlun  were  in  better  agreement  with  the  exact  slope  values 
when  compared  point  by  point,  as  discussed  below. 

33.  Table  2 presents  a comparison  of  the  exact  slope  values  with 
those  calculated  with  the  three  algorithms.  The  comparison  is  made  by 
taking  differences  between  the  exact  and  the  various  algorithmic  values 
on  a grid  point-by-point  basis.  'ITie  values  expressed  in  this  table  rep- 
resent the  percent  probability  that  the  difference  between  the  algo- 
rithmic and  exact  slopes  fell  within  a specific  slope  difference  range. 
Arrays  used  for  all  comparisons  contained  100  by  100  grid  points. 

3^.  Inspection  of  Table  2 indicates  that  the  slopes  computed  with 
the  surface  algorithm  compare  more  favorably  with  the  exact  surface 
slopes  than  do  the  slopes  computed  with  the  other  two  algorithms.  From 
the  point-by-point  differences  computed,  93.8  percent  fell  within  the 
difference  range  of  0 to  <1  deg,  meaning  that  93.8  percent  of  all 
surface  algorithm  values  were  within  1 deg  or  less  of  the  exact  values, 
and  the  remaining  surface  algorithm  values  were  between  1 and  2 deg. 


Nummary  of  Results 


35.  Since  the  primary  objective  of  this  study  was  to  develop  a 
generally  applicable  automated  procedure  to  construct  baseline  slope 
data  (i.e.  slope  maps),  to  support  a variety  of  data  requirements  at 
Army  installations,  no  attempt  was  made  to  evaluate  the  accuracy  with 
which  the  three  algorithms  would  calculate  the  aspect  angles  of  the 
generated  surfaces.  Therefore,  the  only  criteria  for  determining  the 
final  selection  of  an  algorithm  were  the  cost  and  accuracy  of  the  slope 
calculation. 

36.  On  the  basis  of  the  comparisons  made  above  and  the  fact  that 
the  total  computer  computational  cost  did  not  differ  significantly  from 
one  algorithm  to  the  next,  the  surface  algorithm  was  selected  as  the 
"best"  algorithm  for  constructing  baseline  slope  maps  and  was,  there- 
fore, Incorporated  into  the  automated  procedure. 


PART  III:  AN  AUTOMATED  PROCEDURE  FOR  CONSTRUCTING 

SLOPE  MAI’S 


37.  A general  flow  diagram  delineating  the  functional  steps  in 
the  automated  procedure  (computer  program)  for  constructing  slope  maps 
is  presented  in  Figure  11.  In  the  operation  of  the  computer  program, 
named  SLOPEKAP,  the  program  basically  accepts  topographic  data  prepared 
in  advance,  calculates  and  classes  slope  values,  and  then  processes 
these  data  into  optional  forms  of  output.  A computer  listing  of  SLOPE- 
MAT'  and  a glossary  of  program  variables  are  presented  in  'Volume  II  of 
this  report. 


Topographic  Input  Data 

38.  Two  acceptable  forms  of  topographic  data  can  be  input  to 
SLOPJXAP,  elevation  or  vector  data.  The  basic  elevation  data  needed  can 
be  obtained  (see  paragraphs  6 and  T)  from  field  survey  data,  aerial 
photography,  or  topographic  maps.  The  alternative  vector  form  of  in- 
putting topographic  data  is  provided  to  reduce  the  amount  of  computer 
processing  time  during  subsequent  reclassification  runs  on  the  same 
terrain  data.  For  example,  if  a given  terrain  surface  is  to  be  classed 
by  different  slope  class  ranges  (see  paragraphs  U8-50)  several  times, 

it  is  faster  and  cheaper  to  let  SLOPQiAP  generate  a vector  grid  array 
(see  paragraphs  62  eind  63)  during  the  first  run  with  an  elevation  grid 
array  and  then  use  the  vector  grid  array  instead  of  the  elevation  data 
in  all  subsequent  runs.  This  prevents  the  computer  from  having  to 
duplicate  numerous  identical  calculations.  The  general  instructions 
for  preparing  these  two  input  forms  are  describea  below. 

Elevation  data 

39.  Hie  topographic  surface  must  be  represented  by  an  elevation 
grid  array  (Figure  12),  each  element  of  which  is  the  elevation  (in 
metres)  of  the  topographic  surface  at  that  grid  position.  The  two 
planimetric  dimensions  are  deslgi.ated  X and  Y , and  the  elevation 
of  the  surface  Z is  designated  as  the  vertical  dimension  for  the 
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Figure  11  (sheet  2 of  2) 
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Figure  12.  Elevation  grid  array  representing  the 
topographic  surface 

rectilinear  coordinate  system.  (The  input  format  for  the  elevation  grid 
arrays  is  given  in  paragraph  Ul. ) Also,  in  Figure  12,  the  notation 
NYRE  indicates  the  number  of  rows  in  the  matrix  and  NXRE  indicates 
the  number  of  columns.  (See  Volume  II,  Table  II-2. ) A dimension 
statement  restricts  the  grid  array  to  500  columns  or  less  (NXRE  ^ 500). 
Since  X and  Y values  are  xuiderstood  from  the  position  of  the  element 
in  the  array,  only  the  elevation  values  need  be  recorded.  To  provide 
the  user  with  the  capability  of  designating  certain  areas  as  being 
outside  an  area  of  interest,  elevation  values  at  each  grid  point  cor- 
responding to  the  outside  areas  can  be  given  a code  value  of  -999-  This 
will  cause  the  computer  to  assign  values  of  0 to  the  slope  vector  com- 
ponents and  to  insert  a blank  character  into  the  classed  slope  value 
array . 

^40.  The  elevation  grid  array  must  be  prepared  for  a sequential 
computer  file,  such  as  magnetic  tape,  to  be  compatible  with  the  input 
format  of  SLOPEMAP.  The  array  must  be  preceded  on  the  sequential  file 
by  appropriate  identification  records,  which  must  include  the  following: 
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Identification  number.  The  user  must  provide  an 
integer  identification  number  not  to  exceed  10  digits. 

b.  Data  scale.  The  scale  of  the  elevation  data  (DATADC) 
given  in  fractions  of  a metre  (e.g.,  0.01  implies  that 
the  elevation  values  on  the  sequential  file  are  in 
centimetres ) . 

c_.  Grid  size.  The  number  of  rows  (NYRK)  and  columns  (NXRE) 
defining  the  grid  size  in  integer  units. 

d.  Grid  spacing.  The  spacing  (D)  between  the  elements  of 
the  array  is  constant  within  any  one  grid  array,  but  the 
value  must  be  specified  by  the  user.  Values  are  typically 
chosen  in  the  range  of  25  to  100  m. 

^1.  After  the  identification  records  are  placed  on  the  sequential 
file,  the  elevation  grid  array  must  follow  one  row  at  a time,  with  each 
row  treated  as  a single  record.  Thus,  the  final  file  contains  the 
following  records  (the  input  format  for  each  record  is  provided  in  the 
comment  statements  below): 


Record 

1 


2-7 


8 


9 

through 
(NYRE  + 8) 


Comments 

Reserved  for  user's  identification  number. 

FORMAT  (110,  122X) 

Reserved  for  user's  descriptive  information. 

Usually  contains  file  identification  and  the 
location  and  description  of  the  site.  FORMAT  (22A6) 

Contains  the  variable  values  of  DATAGC  , NXRE  , 
NYRE  , and  D , respectively.  FORMAT  (F5.2,  5X,  2l6, 
F6.1) 

Contains  elevation  values  (binary  integers).  The 
number  of  records  from  9 through  (NYRE  + 8)  is 
equal  to  the  number  of  rows  in  the  grid  elevation 
array.  Each  record  contains  the  same  number  of 
elevation  values  as  there  are  columns  in  the  array 


1*2.  Magnetic  tapes  of  digital  elevation  data  can  be  prepared  at 
the  Defense  Mapping  Agency  Topographic  Center  (DMATC)*  or,  in  excep- 
tional instances,  at  the  WES.  An  index  of  currently  available  DMATC 
large-scale  (l:2l*, 000-1: 62, 500)  digital  terrain  elevation  data  on 


* Defense  Mapping  Agency  Topographic  Center,  6500  Brooks  Lane,  Washing- 
ton, D.  C.  20315. 
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military  in:;taliatioris  in  tiie  United  States  is  contained  in  Table  3. 
However,  the  DMATC  data  tapes  are  not  directly  compatible  with  SLOF'H’iAP. 
Therefore,  an  interface  program  (available  from  the  WES)  must  first  be 
used  to  reformat  the  liflATC  tapes  before  these  data  are  acceptable  to 
SLOPIMAP.  The  cost  of  preparation  of  data  tapes  depends  on  the  size  of 
the  area  to  be  digitized,  the  topographic  complexity  of  the  area  of 
interest,  and  the  nature  of  the  primary  data  source  (i.e.,  field  survey 
data,  existing  topographic  maps,  or  aerial  photographs).  Inquiries  con- 
cerning time  and  cost  estimates  can  be  addressed  to  the  DMATC  or  the  WES. 
Vector  data 

U3.  As  stated  in  paragraph  37,  the  representation  of  the  topo- 
graphic surface  can  also  be  input  to  GLOPIMAI'  in  the  form  of  components 
of  the  slope  vectors  (Equation  l).  Instead  of  elevation  values 

representing  each  grid  location,  the  vector  components  f^  - -3f/3x 
and  f^  = -3f/9y  represent  the  topographic  surface  at  each  grid  lo- 
cation in  the  array  (Figure  13). 
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Figure  13.  Vector  grid  array  representing  the  topographic  surface 

I4I4.  A sequential  computer  file  (e.g.  magnetic  tape)  must  be 
prepared  to  make  the  vector  component  data  acceptable  to  SLOPE34AF.  The 
identification  and  vector  grid  array  records  must  be  placed  on  the 

3»* 


1 


- ‘.'i- 


— 


1 


sequential  file  by  exactly  the  srune  procedure  previously  discussed  for 
the  elevation  (-’.rid  array  data  (see  parafi;rapiis  1*0  find  iil).  'Idle  one 
exceptioti  is  that,  instead  of  single  elevation  valuiis,  records  9 thr'iugh 
( NYRE  + 8)  contain  vector  component  pairs  ol'  values  as  binary  integers 
multiplied  by  100,000.  Tiie  number  of  records  from  9 through  (HYRE  ♦ 8) 
again  equals  the  number  of  rows  in  the  grid  array,  but  each  record 
contains  twice  tiie  number  of  values  as  there  are  columns  in  the  array. 

Card  Input  Data 

'fhe  execution  of  SLOPEMAP  is  controlled  by  the  card  input 
data.  These  cards  are  prepared  with  information  pertaining  to  the  form 
of  topographic  input  data,  calculation  units,  slope  class  ranges,  and 
the  desired  output  options.  Each  type  of  card  input  is  discussed  in 
this  section  and  examples  are  shown. 

Form  of  topographic  input  data 

k6.  The  first  step  in  the  execution  sequence  of  SLOPlMAi  (Fig- 
ure 11)  is  to  establish  the  form  of  topographic  input  data  to  be  used, 
i.e.  elevation  grid  array  values  or  vector  component  values.  The  se- 
lection is  made  by  entering  either  the  term  ELEVAT  or  VECTOR  on  the 
first  input  card  (Figure  iha). 

Calculation  units 

hj.  'file  second  input  requirement  is  to  establish  the  units  in 
which  the  slope  values  are  to  be  calculated,  i.e.,  slope  values  expressed 
in  degrees  or  tangent  values  of  the  angle.  All  subsequent  input  and 
output  data  are  then  in  the  chosen  units.  The  selection  is  made  by 
entering  into  the  second  card  either  the  term  DECREE  or  TANGEN 
(Figure  ll*b). 

Slope  class  ranges 

1*8.  The  slope  class  ranges  are  read  into  the  program  on  succeed- 
ing input  cards.  'The  parameters  punched  into  card  3 are  the  lower  and 
upper  bounds  of  the  first  slope  class  range.  Each  succeeding  pair  of 
values  encountered  on  cards  1*,  5,  6,  etc.,  define  the  second,  third, 
fourth,  etc.,  slope  class  ranges  (Figure  lUc).  The  third  card  should 
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OUTPUT  OPTIONS 


contuin  two  neRative  vaiueo  separated  by  commas  if  no  class  values  are 
to  be  output . 

149.  The  progrim  format  interprets  all  input  data  cards  for  slope 
class  range  by  a "list-directed"  format  statement  (statement  1 of  the 
main  program  in  Table  ll-l.  Volume  II),  The  list-directed  form  of  en- 
tering data  implies  that  the  number  and  types  of  list  variables  directly 
determine  the  way  in  which  the  comi)uter  assigns  values  to  the  variables. 
This  offers  a much  more  flexible  form  of  inputting  the  slope  class 
ranges  as  compared  to  a fixed  columnar  form  of  format  statement.  The 
one  requirement  of  the  list-directed  form  of  inputting  is  that  the 
variable  values  on  each  card  be  separated  from  each  other  by  at  least 
one  blank  column  or  a comma. 

50.  Note  that  the  last  class  range  card  in  Figure  lUc  has  a pair 
of  negative  values.  Any  negative  number  will  end  the  reading  of  classed 
values  and  cause  control  to  be  passed  on  to  the  class  range  editor.  The 
editor  checks  for  possible  errors  and  outputs  a message  to  the  printer 
if  errors  are  detected.  The  last  card  containing  negative  values  is 
still  required  even  if  no  classed  values  are  to  be  output. 

Desired  output  options 

51.  The  last  data  card  (Figure  lUd)  read  by  the  program  contains 
the  output  control  options.  Table  It  contains  a list  of  all  the  control 
values  and  short  descriptions  of  their  functions.  Each  output  form  is 
explained  more  fully  in  the  following  section. 

Output  of  ELOPEMAl^ 

52.  SLOPEMAP  can  be  directed  by  the  output  opt.' on  card  to  process 
the  calculated  and  classed  slope  information  in  one  or  more  of  four 
modes.  This  section  discusses  the  four  output  modes  available  as  well 
as  two  additional  processing  options  of  the  calculated  and  classed  slope 
data. 

Output  modes 

53.  Swath  dump.  The  rows  of  the  calculated  and  classed  value 
arrays  are,  in  general,  too  long  to  be  output  on  one  line  of  the  printer. 
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The  jiro^ram  solves  this  by  calling  the  subroutines  CALilWH  or  CLilGWH 
to  print  vertical  swaths  of  values.  Each  swath  corresponds  to  a swath 
taJ<en  from  the  calculated  or  classeci  arrays  and  is  restricted  to  a 
maximum  of  132  characters  per  line.  After  the  first  swath  has  been 
written  by  the  printer,  the  program  begins  the  next  swath  and  continues 
until  the  entire  set  of  array  elements  has  been  printed.  For  display 
purposes,  the  separate  swaths  can  be  laid  side  by  side  and  assembled 
manually  to  form  a mosaic  slope  map  of  computer  calculated  or  classed 
slope  values  as  shown  in  Figure  I5. 

5^^.  Magnetic  tapes.  The  user  has  the  option  of  storing  calcu- 
lated or  classed  slope  values  onto  magnetic  tape  in  binary  integer  form. 
The  program  calls  subroutine  TaPCAL  to  output  a tape  of  calculated 
values  or  subroutine  TAPCLS  to  output  a tape  of  classed  values.  ITie 
calculated  slope  values  are  recorded  on  tape  following  the  descriptive 
information  copied  from  the  input  topographic  data  tape.  Classed  values 
are  recorded  following  the  descriptive  information,  the  number  of  classes 
(HCLAS),  and  the  previously  selected  class  ranges.  Because  of  inherent 
round-off  errors  during  integer  arithmetic  performed  ‘by  the  computer, 
SLOPEMAP  was  programmed  to  calculate  slope  values  at  a scale  100  times 
greater  than  actual  values  for  degree  units  and  100,000  times  greater 
than  actual  values  for  tangent  units.  Therefore,  the  calculated  values 
output  to  magnetic  tape  represent  slope  values  multiplied  by  these 
amounts.  The  tapes  contain  the  following  records  (the  output  format  for 
each  record  is  provided  in  the  comment  statement  below): 
li.  Calculated  slope  values. 

Record  Comments 

1 Contains  the  identification  number  copied 

from  the  input  topographic  data  file. 

FORMAT  (I10,122X) 

2-7  Contains  the  descriptive  information  copied 

from  the  input  topographic  data  file. 

FORMAT  (22A6) 

8 Contains  variable  values  DATASC  , NXRE  , 

NYRE  , and  D , respectively.  FOIFIAT 
(F5.2,5X,2I6,F6.1) 
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Figiire  15.  An  assembled  swath  dump  of  classed  values  for  a generated  cuspf unction  surface 
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Record 


Comments 


9 Contains  calculated  slope  values  (binary 

through  integers).  The  number  of  records  from  9 

(NYRE  + 8)  through  (NYRE  + 8)  is  equal  to  the  number 

of  rows  in  the  original  input  topographic 
grid  array.  Each  record  contains  the  same 
number  of  slope  values  as  there  were 
columjis  in  the  input  array.  Each  slope 
value  has  been  multiplied  by  either  100 
(for  degrees)  or  100,000  (for  tangent  of 
of  the  angle ) 


b.  Classed  slope  values. 

Record Comments 


1 Contains  the  identification  number  copied 

from  the  input  topographic  data  file. 
FORMAT  (I10,122X) 

2-7  Contains  the  descriptive  information 

copied  from  the  input  topographic  data 
file.  FORMAT  (22A6) 

8 Contains  variable  values  of  DATASC  , 

NXRE  , NYRE  , and  D , respectively. 
FORMAT  (F5.2,5X,2I6,F6.1) 

9 Contains  the  number  of  slope  class  ranges 
NCLAS  (in  binary) 


10 

through 
(IJCLAG  + 9) 


Contains  the  lower  and  upper  bounds  of 
each  class  range.  The  number  of  records 
from  10  throiigh  (NCLAo  + 9)  is  equal  to 
the  number  of  slope  class  ranges.  Each 
record  contains  a pair  of  numbers  (binary) 


(NCLAS  + 10) 
through 
(NYRE  + NCLAS 
9) 


Contains  classed  slope  values  (binary  in- 
tegers). The  number  of  records  (NCLAS  + 10) 
through  (NYRE  + NCLAS  +9)  is  equal  to  the 
number  of  rows  in  the  original  input  topo- 
graphic grid  array.  Each  record  contains 
the  same  number  of  classed  values  as  there 
are  columns  in  the  input  array 


55-  Punch  cards.  Calculated  or  classed  values  can  be  output  into 


punched  cards  by  the  subroutines  CRDCAL  and  CRDCLS  , respectively. 
Descriptive  information  contained  on  the  elevation  data  tape  is  first 
punched  by  the  selected  subroutine.  The  calculated  values  follow  the 
descriptive  information  in  the  case  of  CRDCAL  ; but  in  the  case  of 
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CRDCLL)  , the  number  of  clacses  and  class  ranges  rre  punched  prior  to  the 
classed  values  being  output  in  integer  form.  The  card  decks  contain  the 
following  information  (the  output  format  for  each  card  is  provided  in 
the  comment  statements  below): 

a.  Calculated  slope  values. 

Card Comments 


1 Contains  the  identification  number  copied 

from  the  topographic  data  file  and  a card 
sequence  number.  FORMAT  ( 110 ,60X , 110 ) 

2-7  Contains  the  descriptive  information 

copied  from  the  input  topographic  data  file 
and  a card  sequence  number.  FORMAT  (11A6, 

Ux,iio) 


8 Contains  variable  values  DATAOC  , NXRE  , 

NYRE  , D , and  a curd  sequence  number. 
FORMAT  (F5.2,5X,?l6,F6.1,i42X,I10) 


9 

through 

NYRE(WXRK/7) 

+ 8 or  NYRE 
(NXRE/T  + 1) 

+ 8 


Contains  calculated  slope  values  and  a 
card  sequence  number.  The  maximum  number 
of  calculated  slope  values  per  card  is  7- 
If  (NXRE/7)  is  an  integer  value,  then  the 
total  number  of  cards  containing  calcu- 
lated values  equals  NYRE(NXRE/7)  • How- 
ever, if  (NXRE/7)  is  not  an  integer  value, 
then  the  number  of  cards  containing  cal- 
culated value  equals  NYRE(NXRE/7  + 1 ) . 
Each  output  slope  value  is  either  100 
(for  degrees)  or  100,000  (for  tangent  of 
the  angle)  times  greater  than  their  actual 
magnitudes.  P’ORMAT(10I7,I10) 


b.  Classed  slope  values. 


Card  Comments 

1 Contains  the  identification  number  copied 

from  the  input  topographic  data  file  and  a 
card  sequence  number.  FORMAT  (I10,60X,I10) 

2-7  Contains  the  descriptive  information  copied 

from  the  input  topographic  data  file  and  a 
card  sequence  number.  FORMAT  (11A6,UX,I10) 

8 Contains  variable  values  UATASC  , NXRE  , 

NYRE  , D , and  a card  sequence  number. 
FORMAT  (F5.2,5X,2l6,F6.1,l42X,I10) 
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Comments 
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9 Contains  the  number  of  slope  class  ran^'es 

(NCLAC)  and  a card  sequence  number.  FORMAT 
(I10,60X,I10) 

10  Contains  lower  and  upper  bounds  of  each 

through  class  rfinge  and  a card  sequence  n'jmber. 

(NCLAI3  9)  The  number  of  cards  from  10  through  (NCLAi' 

+9)  is  equal  to  the  number  of  slope  class 
ranges.  Each  card  contains  three  numbers. 
FORMAT  (2F10.3,50X,I10) 

(iiCLAll  + 10)  Contains  classed  slope  values  and  a card 

througii  sequence  number.  Tlie  maximum  number  of 

NYRE(NXRE/35)  classed  slope  values  per  card  is  35.  If 

+ (NCLAO  + 9)  (liXRE/35)  is  an  integer,  then  the  total  num- 

or  ber  of  cards  containing  classed  values  equals 

NYRE(NXRE/35  NYRE(:iXRE/35 ) • However,  if  (NXRE/35)  is  not 
+ l)  + (NCLAO  an  integer,  then  the  number  of  cards  containing 
+ 9)  classed  values  equals  WYRE( NXRE/35  1 ) . 

FORMAT  (3512,110) 

56.  Drum  and  CRT  plots . If  the  user  wants  a graphic  display  of 
the  boundaries  between  delineated  classed  values,  he  must  choose  the 
option  to  output  either  a drum  plot  or  a CRT  plot.  Each  outputs  a slope 
map  with  lines  separating  different  classed  values.  Figures  8-10  show 
example  drum  plots  of  slope  class  boundaries.  Descriptive  information 
contained  on  the  input  topographic  data  tape  is  also  plotted. 

57.  The  scale  for  outputting  the  drum  plot  must  be  set  so  that 
the  X and  Y dimensions  of  the  plot,  XLEN  and  YLEN  , do  not  ex- 
ceed 6U  and  2^  in.  (limitation  of  WES  drum  plotter),  respectively. 

Therefore,  the  user  must  carefully  select  the  drum  plot  scale  value, 

PLSCL  , before  running  the  program.  The  following  set  of  equations  can 
be  used  to  check  selected  plot  scales  to  ensure  that  XLFT'J  and  YIJTi 
do  not  exceed  their  limits; 

XLEN  (inches)  = (PLSCL)x(  39.  37007  )>‘(D)x(NXRE-l ) 

Y1J:n  (inches)  = (PLSCL)>‘(39. 37007)x(D)x(NYRE-l) 

where 

D = elevation  grid  spacing  in  metres 
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NXHE 

nyre; 


= number  of  grid  points  along  the  X and  Y ajces, 
respectively 

The  CRT  plot  does  not  require  a scale  input  value  but  generates  its 
own  at  run  time.  The  program  selects  and  outputs  a scale  that  will 
allow  the  CRT  plot  to  be  output  as  large  as  possible,  but  within  its 
maximum  size  of  IT  by  11  in.  (limitation  of  WES  CRT  plotter).  This 
feature  makes  available  a reliable,  inexpensive,  and  juick  graphic  aid 
for  use  during  the  early  stages  of  slope  map  construction.  After  out- 
put needs  have  been  established,  the  more  costly,  but  more  accurate, 
drum  plot  option  can  be  used. 

Additional  processing  options 

58.  Smoothing.  SLOPEMAP  offers  the  capability  of  further 
processing  the  calculated  values  and  classed  values  in  two  additional 
ways  before  output.  The  decision  to  use  one  or  both  ways  is  made  on 
the  output  option  card  in  card  field  5 (Table  U and  E’igure  lUd). 

59"  If  card  field  5 contains  a 0,  no  further  processing  will  be 
perfomed  on  either  the  calculated  values  or  classed  values.  If  the 
fifth  position  contains  a 1,  the  subroutine  SMOOTH  is  called  after 
the  program  finishes  calculating  slope  values.  This  routine  calculates 
an  average  X component  and  Y component  from  the  slope  vectors  sur- 
rounding a given  grid  point.  Included  in  the  average  are  the  nearest- 
neighbor  and  the  next -nearest-neighbor  slope  vector  components.  The 
slope  components  of  the  grid  point  itself  are  excluded.  The  average 
components  are  then  used  to  calculate  a new  slope  value  for  the  given 
grid  point,  and  this  new  value  permsuiently  replaces  the  old  one.  The 
process  has  the  effect  of  smoothing  or  averaging  out  rapid  slope  changes. 

60.  If  the  fifth  column  contains  a 2,  the  subroutine  ISLAND  is 
called  after  the  calculated  slope  values  have  been  classed.  Tliis  rou- 
tine removes  isolated  classed  values  from  the  classed  value  array.  For 
a given  grid  point,  this  subroutine  examines  the  nearest  and  next- 
nearest  neighbors  to  determine  if  the  given  point  is  isolated,  i.e.  no 
equal  class  contiguous  to  it.  If  the  isolated  value  does  not  differ  by 
more  than  one  class  value  from  its  neighbors,  it  is  replaced  by  the 
closest  class  value.  If  the  isolated  value  falls  equally  between  a 
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higher  (iiid  a lower  neighbor  value,  it  takes  the  lower  value,  ITie  sub- 
routine does  not  alter  the  isolated  value  if  its  neighbors  differ  by  two 
or  more  classed  values.  The  removal  of  isolated  islands  existing  in  the 
classed  value  array  can  have  the  effect  of  removing  features  that  would 
result  in  poor  delineation  of  the  slope  regime  of  the  topographic  surface. 

61 . Both  of  the  two  options  can  be  selected  to  operate  on  the 
data  sequentially  by  choosing  a value  of  3 I'or  the  fifth  card  field  of 
the  output  option  card.  The  subroutine  SMOOTH  will  be  executed  prior 
to  subroutine  ISl^ND  . 

62.  Vector  components.  The  slope  algorithm  makes  use  of  the  X 

and  Y components  of  the  slope  vector  to  calculate  slope  values.  If 
directional  information  is  needed  for  later  use  (e.g.  as  a subsequent 
input  tape  to  SLOPEMAl’),  then  the  vector  components  can  be  output  onto 
magnetic  tape.  Since  the  slope  vector  components  are  calculated  and 
stored  within  the  computer  at  100,000  times  their  actual  magnitudes  as 
was  similarly  done  for  the  slope  values  (see  paragraph  the  com- 

ponent values  are  output  to  magnetic  tape  at  this  scale. 

63.  The  three-dimensional  vector  fJ  can  be  reconstructed  to  rep- 
resent the  slope  vector  (or  vector  normal  to  the  surface)  at  any  given 
grid  location  (see  paragraph  lb).  The  slope  vector  is  then  written 

N = (XCOMP/100000)  i + ( YCOMP/100000 ) J + k 

where  XCOMP  = 100000( -3f /3x ) and  YCOMP  = 100000(-3f/9y ) are  the 
vector  components  stored  on  the  slope  vector  output  tape.  Tlie  tape 
contains  the  following  records  (the  output  format  for  each  record  is 
provided  in  the  comment  statements  below); 

Record  Comments 

1 Contains  the  identification  number  copied  from  the 

input  topographic  data  file.  FORMAT  (I10,122X) 

2-7  Contains  the  descriptive  information  copied  from 

the  input  topographic  data  file.  FORMAT  (22A6) 

8 Contains  variable  values  DATASC  , NXRE  , NYRE  , 

and  D , respectively.  FORMAT  (F5.2,5X,2I6,F6.1) 


Hecord 


Comment r 


9 

throiiKh 
(HYRE  + 8) 


Contains  X and  Y sloije  vector  component  pairs 
(binary  integers).  Each  component  has  been  multi- 
plied by  .100,000.  The  number  oT  records  from  9 
through  (NYRE  +8)  is  equal  to  the  number  of  rows 
in  the  grid  array,  and  each  record  contains  twice 
the  number  of  values  as  there  are  coliomns  in  the 
array 


Error  Messages 


6U.  Error  messages  are  incorporated  into  the  computer  program  to 
assist  the  user  in  identifying  input  errors.  These  error  messages  are 
listed  below;  comment  is  provided  if  necessary. 

a.  "The  input  topographic  data  were  not  specified  in  an 
acceptable  foi’m.  The  only  choices  available  to  the  user 
are  ELEVAT  or  VECTOR  . This  error  causes  immediate 
program  termination.  Please  correct  the  input  before 
resubmitting. " 

b.  "The  input  calculation  units  have  not  been  properly 
selected.  Only  choices  are  DEGREE  or  TANGEN  . This 
error  causes  iiraiiediate  program  termination.  Please 
correct  the  input  before  resubmitting." 

£.  "The  number  of  input  slope  classes  exceeds  the  present 
allowed  maximum  36.  If  this  is  not  an  error  and  you 
think  the  output  graphics  will  still  look  reasonable  with 
so  many  classes,  you  will  have  to  change  the  dimension  of 
the  variable  CLASS  and  the  value  of  the  constant 
MAXCLASS  and  resubmit  the  run." 

"The  class  ranges  on  the  input  cards  should  not  overlap. 
This  error  causes  immediate  program  termination.  Please 
correct  the  input  before  resubmitting." 

e_.  "There  are  duplicate  cards  in  the  slope  class  range  input 
cards.  A missing  card  is  suspected.  This  error  causes 
immediate  program  termination.  Please  correct  the  input 
before  resubmitting." 

£.  "The  class  range  input  cards  are  out  of  sequence.  This 
error  causes  immediate  program  termination.  Please 
correct  the  input  before  resubmitting." 

"To  output  a magnetic  tape  with  vector  components,  the 
control  value  on  the  output  option  card  must  be  0 or  1." 

h.  "The  control  valuas  appearing  in  the  first  three  card 
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rie.ldis  on  the  output  option  ear<l  must  bo  in  the  ranp.o  0, 

J,  2,  or  3.  Please  correct  the  output  option  data  card." 

i_.  "The  control  value  for  the  plot  option  must  be  0,  1,  or 
2.  Please  correct  the  output  option  data  card." 

j_.  "The  control  value  for  the  smoothing  options  on  the 

output  option  card  must  be  0,  i , 2,  or  3-  Please  correct 
the  output  option  data  card." 

It.  "No  outp'ut  options  have  been  selected.  Please  correct 
the  output  option  data  card.  " 

_1_.  "You  asked  for  an  output  of  classed  values  but  did  not 

define  class  rariges.  This  error  causes  immediate  program 
termination.  Please  correct  the  input  before  resubmitting. 

m.  "You  are  requesting  a drum  plot  with  dimensions  greater 
than  60  » 24  in.  (X  by  Y dimensions).  Please  recheck 
your  plot  option  scale." 

t£.  "The  grid  array  of  elevations  has  an  insufficient  number 
of  data  points  to  perform  a slope  calculation;  the  array 
size  must  be  at  least  a 4 x 4 matrix." 


Job  Control  Language  for  SLDPEMAP 


65.  A complete  set  of  Job  control  language  (JCL)  cards  and  an 
exfimple  set  of  data  necessary  to  run  SLOPEMAi’  on  the  WES  Honeywell 
CARDIN  time— sharing  (T/S)  system  are  presented  in  Figure  I6.  A de- 
tailed description  of  each  card  function  is  provided  in  the  Honeywell 
reference  manuals. The  Honeywell  G-635  T/S  system  (system  release  8 
update  3)  was  used  throughout  the  development  of  the  program.  Thus, 
the  JCL  cards  are  necessarily  tailored  to  this  system;  however,  they 
should  still  provide  guidance  in  setting  up  JCL  cards  or  their  equiva- 
lent for  other  computer  system.s  such  as  those  that  might  be  available 
on  military  installations. 

66.  It  should  be  noted  that  a typical  run  does  not  require  the 

use  of  all  cards.  For  example,  only  the  $: USE: DRUM  or  the  $: USE: CRT 
should  be  inserted  if  a drum  or  CRT  plot  is  requested  by  the  user.  The 
same  is  true  of  any  of  the  $:TAPET:  cards  and  the  message  cards 
$:MSG2:  following  them.  Only  if  output  is  requested  on  these  devices 

should  these  cards  be  inserted.  Also,  only  one  topographic  input  data 
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Figure  l6.  Complete  set  of  JCL  cards  to  run  SLOPEMAP  from  a source 
deck  on  the  Honeywell  G-635  system 
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device  is  required.  Tlierel’ore,  the  unused  input  file  cards  (e.g. 

and  ^rPRMEL:  ) must  be  deleted.  Correspondingly,  in  the 
control  data  cards,  only  one  form  of  input  topographic  data  and  one  tyj>e 
of  calculation  units  should  be  chosen.  A list  of  file  codes  and  their 
respect ! /e  input  and  output  file  descriptions  are  presented  in  Table  5 
to  help  the  user  assemble  a JCL  deck. 

67.  If  a non-WKo  user  of  SLOPEMAP  cannot  provide  the  function 
implied  by  the  $:SELECT:  card  (Figure  I6),  he  should  replace  this 
card  with  the  SLOPEMAP  source  statement  cards  as  listed  in  Table  II-l 


of  Volume  II.  For  a WES  user,  the  SLOPEMAP  program  should  reside  as  a 

binary  coded  decimal  (BCD)  file  under  his  user  identification  number. 

In  both  cases,  execution  follows  compilation  of  the  source  statements. 

68.  To  avoid  the  necessity  of  compilation  before  every  computer 

run,  the  WES  user  can  take  advantage  of  an  existing  magnetic  tape  that 

contains  a compiled  and  loadable  version  of  SLOPEMAP.  This  tape, 

called  a system  loadable  product  tape,  was  produced  by  the  Honeywell 
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Production  Library  Generator  ’ ’ and  is  maintained  for  WES  users.  The 
JCL  cards  to  run  this  version  of  SLOPEIiAP  are  presented  in  Figure  17. 
The  same  restrictions  with  respect  to  JCL  cards  still  apply  with  the 
exception  of  the  $;PRODUCT;  card.  Should  the  WES  user  choose  to 
output  either  a CRT  plot  or  a drum  plot,  he  must  use  only  a $: PRODUCT: 


SLOPEC  or  a $; PRODUCT : SLOPED  capd.  A current  tape  number  correspond- 


ing  to  the  system  loadable  product  tape  can  be  obtained  by  addressing 


an  inquiry  to  the  WES. 

69.  The  SLOPEMAP  program  is  available  to  any  authorized  user  with 
a remote  terminal,  such  as  personnel  at  military  installations.  To  ob- 
tain authorization,  a user  must  first  obtain  a WES  computer  account 


number  and  information  on  accessing  the  system.  Authorization  for  use 
of  the  SLOPEMAP  program  on  the  WES  computer  can  be  obtained  by  contact- 
ing the  WES.  Telephone  inquiries  should  be  made  to  the  author  at  com- 
mercial telephone  601-636-3111,  or  AUTOVON  1*35-1680,  extension  WES,  or 
FTS  51*2-3111. 
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aas  J I DENT:  USD»I  CW  STRir^E 
30$:t»POnJCT:  SLOPED 
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A0S:omfL:  l0#R#L#USERID/nATAC#R 
90$:TAPEr7:Il,XICW#(TAPE  NUMBER) 

1 00  S:  TAPE?;  20#  X2DR#  # # # CALCULATED-  WLUES 

II0$:MSG2;SA'/E  20#  STRUVE# USERID.  CM-CILATED-VALUES  OF  E-OreS 
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Figure  17.  Complete  set  of  JCL  cards  to  run  SLOPEMAP  from  a 
system  loadable  product  tape  on  the  Honeywell  G-6l5  system 
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PART  IV:  DEMONilTRATION  OF  THE  AUTOMATED  PROCEDURE 
FOR  CONSTRUCTING  SLOPE  MAPS 


70.  To  demonstrate  how  baseline  slope  data  can  be  obtained  usin>^ 
the  program  SLOPEMAF,  two  sites  were  chosen,  one  located  in  the  Middle 
East  (Israel)  and  the  other  was  located  at  an  Army  installation  (Hunter 
Liggett  Reservation,  California).  These  sites  were  selected  because 
they  contained  two  different  contour  intervals,  20-m  (65.6-ft)  contours 
for  the  Middle  East  site  and  UO-ft  contours  for  the  Hunter-Liggett  site 
The  map  scales  of  the  two  sites  were  1:50,000  for  the  Middle  East  site 
and  1:25,000  for  the  Hunter-Liggett  site.  The  approach  taken  was  to 
choose  example  outputs  from  the  available  optional  outputs  and  then  to 
evaluate  them  with  respect  to  manually  produced  slope  maps.  The  pro- 
cedures for  producing  slope  maps  of  both  sites  are  discussed  and  evalu- 
ated in  the  following  paragraphs. 

Middle  East  Site  (l;  50,000  -Map  Scale) 

71.  The  procedural  steps  followed  to  produce  slope  maps  of  the 
Middle  East  site  were  to: 

a.  Construct  input  elevation  grid  array. 

Edit  and  examine  elevation  grid  array. 

£.  Select  forms  of  desired  output. 

d.  Prepare  input  data  and  JCL  cards  necessary  to  execute 
SLOPEMAP. 

£.  Execute  SLOPEMAP. 

£.  Construct  slope  maps  from  selected  output  forms. 

Site  elevation  data 

72.  A topographic  map  of  the  selected  Middle  East  site  is  shown 

in  Figure  18.  The  northwest  corner  of  the  site  is  6 km  southeast  of 

Nazareth,  Israel.  The  contour  interval  is  20  m with  an  occasional 

supplementary  10-m  contour  shown  in  the  flatter  areas.  Each  contour 

line  within  the  10-  by  7-km  site  was  digitized  and  processed  by  a WES 

3 

computer  program  that  generates  an  elevation  grid  array.  The  spacing 
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Topographic  map  of  the  Middle  East  site  (map  scale  1:50,000) 
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between  t'l'id  points  was  set  at  a constant  50  m.  'Fhe  size  of  the  sub- 
sequent grid  array  was  201  by  lUi  (NXHE  by  NYRE).  TJie  elevation  values 
themselves  were  stored  as  integer  multiples  (metres). 

73.  Once  the  grid  elevation  array  was  constructed,  it  was  edited 
!Uid  exfunined  for  errors.  A method  frequently  used  at  the  WEO  to  quickly 
detect  gross  misrepresentations  of  the  terrain  surface  is  to  portray  the 
data  by  three-dimensional  perspective  plots.  These  plots  allow  rapid 
visual  examination  of  large  elevation  data  arrays.  The  computer  gene- 
rated perspective  plots  of  the  Middle  East  elevation  grid  array  are 
shown  in  Figure  19.  One  of  the  prominent  features  in  these  plots  is  the 
terracing  effect. in  the  terrain  surface  along  and  parallel  to  contour 
lines.  These  terraces  are  not  actual  terrain  conditions  but  are  pro- 
duced artificially  by  the  elevation  interpolation  algorithm  used  in 
the  computer  program  that  generates  the  grid  data  from  the  contour  data. 

This  effect  is  a function  of  the  grid  spacing  constant  and  diminishes  as 

9 

the  grid  spacing  constant  decreases.  With  a grid  spacing  of  50  m,  how- 
ever, the  perspective  plots  do  show  that  the  Middle  East  elevation  grid 
array  does  approximate  the  terrain  surface  and,  for  purposes  of  providing 
input  topographic  data  to  SLOPEMAP,  is  adequate.  If  the  terracing 
effect  had  been  found  to  invalidate  the  data,  the  grid  spacing  constant 
could  have  been  decreased  and  the  grid  elevation  array  regenerated. 

A better  solution  to  the  terracing  artifacts  would  have  been  to  eliminate 
the  effect  altogether  by  modifying  the  original  interpolation  algorithm. 
The  latter  solution  would  have  required  an  extensive  programming  and 
testing  effort  and  thus  was  considered  out  of  the  scope  of  work  reported 
herein. 

Selection  of  output  forms 

'fh.  In  general,  the  slope  classes  selected  by  the  user  will,  of 
course,  be  problem  dependent.  However,  in  lieu  of  a priori  knowledge  of 
any  class  ranges,  the  user  could  apply  the  same  technique  to  establish  a 
tentative  set  of  slope  class  reinges  as  was  used  in  this  study  to  estab- 
lish a demonstration  set  of  slope  class  ranges.  The  technique  was  to 
run  SLOPiMAi'  with  the  slope  class  ranges  set  at  1-deg  increments  ranging 
from  0 throxigh  90  deg.  A graph  was  then  produced  for  both  the  Middle 
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Figure  19  (sheet  3 of  ^4 ) 


East  site  and  the  Hunter-Liggett  site  (see  jiarugraph  8i)  r.howing  the 
accumulative  percentage  of  the  areas  of  the  sites  as  a function  of  in- 
creasing slope  values  (F'igure  20).  Inspection  of  these  graphs  revealed 
that  the  distribution  of  slopes  were  very  similar.  Therefore,  only  one 
set  of  slope  class  ranges  was  selected  for  the  two  sites  to  demonstrate 
the  slope  classification  option  of  SLOPENLM'.  The  class  ranges  selected 
were  as  follows: 

Class  Hange , deg 

A 0 to  <5 

B 5 to  <10 

C 10  to  <30 

D 30  to  <1*5 

E 1*5  to  90 

75-  The  following  forms  of  output  were  then  chosen  to  display  the 
slope  information  for  this  example: 

&.  Swath  dumps  of  calculated  and  classed  slope  values. 

b.  Magnetic  tape  of  classed  values  only. 

£_.  Punched  cards  with  calculated  values  only. 

CRT  plot . 

e_.  No  smoothing  options. 

£.  Slope  vector  components  on  magnetic  tape. 

The  control  values  necessary  for  execution  of  each  of  these  selected 
forms  of  output  are  obtained  from  Table  1*.  F'or  this  exfimple,  the 
control  values  and  their  proper  order  of  placement  on  the  output  option 
card  (Figure  lUd)  are: 


3 , 2 , 1 , 2 , 0 , 1 , 0.0 

The  drum  plot  scale  value  of  0.0  is  used  only  as  a dummy  value,  since 
the  program  does  not  use  this  value  when  called  to  output  a CRT  plot. 

It  would  have  been  sufficient  to  end  the  string  of  numbers  with  a comma, 
in  which  case  the  program  would  have  substituted  a value  of  zero.  The 
minimum  requirement,  however,  is  at  least  a comma  to  define  the  necessary 
card  fields.  This  is  also  true  if  one  wishes  to  null  (l.e.  zero)  any  of 
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Accumulative  percentages  of  the  areas  of  the  demonstration  sites  as  a function  of  slope 


the  oi:tioii  fields.  For  exfimple,  in  the  present  illustration,  the  output 
option  card  could  have  been  set  up  in  the  following  way: 


3 , 2 , 1 , 2 , , 1 , 


T6.  The  JCL,  data,  and  option  cards  necessary  to  execute  the  ex- 
extunple  with  a system  loadable  product  tape  are  shown  in  Figure  21.  The 
user  should  refer  to  Figure  16  if  he  wishes  to  use  the  $:SELECT:  or 
source  deck  form  of  program  input. 
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20$!  1 E«NT: USERl  IV  SmWE 
30 Si PRODUCT: SLOPBC 
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60S:mPE7tl0#Xia#<'mPE  NIM9ER> 

70$:  TAPETi  21  *X3ER*  GLASSED- ’JflLUES 

80$:MSG2:  SAVE  21/ STRUVE» USERID* OASSED-VflLUES  OF  S-OPES 
90$:  TAPE?:  22/X4DR*  / / * S-OPE-VBCTOR 

100$:MSG2:SAVE  22,  STRUVE*  USERIIV  3-OPE- VECTOR  CCMPONEJITS 

1 10$:  TAPE?:  24,X6DR* , * , CRT-FLOT 

120$:MSG2iSAVE  24*  STRUVE* USERIIV  CRT- PLOT 

130$inLE:30,X?R,  10R 

140$inLE:  31,X8R*10R 

1 50$:  HLE:  32,X9R,  20R 
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190#0, 5 

200#5*10 
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260$:ENDJOB 

Figure  21.  Job  control  language,  data,  and  option  cards  used  to 
produce  the  slope  map  in  Figure  2U  of  the  Middle  East  site 

Output  forms  and  reports 

77.  The  initial  output  reported  on  the  high-speed  printer,  after 
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Job  execution,  gives  a listing  of  the  input  instructions,  selected  cal- 
culation units,  slope  class  ranges,  output  options,  and  descriptive  in- 
formation contained  on  the  input  topographic  data  tape.  The  output 
report  for  the  Middle  East  example  is  presented  in  Figure  22,  If  there 
had  been  an  error  detected  in  the  input  instructions,  one  of  the  error 
messages  listed  in  paragraph  61  would  have  been  printed  in  the  initial 
output  report,  and  the  program  would  have  then  been  terminated. 

T8.  Following  the  initial  output  report  of  input  instructions  are 
the  swath  dump  reports  and  values.  Both  the  calculated  and  classed 
slope  value  reports  contain  the  descriptive  information  stored  on  the 
input  topographic  data  tape.  However,  in  the  classed  slope  value  report 
(Figure  23),  an  additional  listing  containing  a legend  of  the  symbols 
used  for  the  various  class  ranges  is  output.  Also  printed  are  the  num- 
ber of  grid  locations  classed  in  each  class  range  under  the  heading 
COUNT  (e.g.,  class  range  0 to  <5  deg  classed  ll*,U90  grid  locations  as  A). 
Immediately  following  the  descriptive  reports  are  the  swath  dumps  them- 
selves. To  produce  the  swath  dump  slope  maps,  it  is  necessary  to 
assemble  the  Middle  East  swath  dumps  of  calculated  and  classed  slope 
values  (see  paragraph  53).  Once  these  report  and  swath  dump  options 
have  been  satisfied,  no  further  use  is  made  of  the  high-speed  printer. 

The  remaining  output  options  (i.e.  magnetic  tapes,  punch  cards,  CRT 
plots,  and  vector  components)  are  output  according  to  their  respective 
formats  and  output  devices  (see  paragraphs  5^-57  and  62). 

79«  Other  forms  of  analysis  and  display,  which  require  SLOPEMAP 
products  as  input  data,  can  be  developed.  An  example  of  a postprocess- 
ing procedure  used  extensively  in  this  study  was  the  use  of  the  classed 
value  tape  to  produce  slope  map  pictures.  Through  this  method  of 
producing  slope  maps,  the  best  visual  comparisons  between  manual  and 
computer  procedures  can  be  made.  This  picture  process  is  discussed  in 
greater  detail  in  paragraphs  80-82. 

Slope  map  pictures 

80,  The  contents  of  the  magnetic  tape  output  of  classed  slope 
values  (see  paragraph  5*+)  consist  of  binary  integers  that  represent  the 
class  ranges  specified  by  the  input  instructions.  For  example,  integers 

t 
1 

\ 

( 


6l 


Fig’ore  22.  Initial  computer  output  report  giving  input  instructions  and  descriptiv' 
information  on  the  input  topographic  data  tape  for  the  Middle  East  example 


]-5  'ire  stored  oii  the  getioruted  in  the  Middle  Kar.t  example  since 

five  classes  were  defined  on  tlie  dntn  cards.  KqcIi  intef'er  corresponds 
to  a legend  symbol  listed  in  the  clause.!  value  swatii  dump  report  (e.g, 
1-A,  2iB,  ,..5iK). 

81 . After  the  classed  slope  value  tape  was  generated  by  SLOi’lWAi", 
it  was  used  as  an  input  tape  for  producing  slope  maps  in  a picture  format 
on  the  we;:  film  reader/writer.  Figure  2U  is  a slope  map  of  the  Middle 
East  site  generated  on  tills  device.  Note  that  the  terracing  effects  de- 
tected previously  in  the  persj-ective  plots  of  the  elevation  data  (Fig- 
ure 19)  are  discernible  here  in  slope  class  rijinge  B . Also,  since  the 
options  to  smooth  and  remove  isolated  classed  islands  were  not  selected 
output  options,  numerous  isolated  blocks  or  squares  are  observed  that 
tend  to  complicate  portions  of  the  slope  map.  This  may  or  may  not  be  of 
importance,  depending  on  the  type  of  problem  for  which  data  are  required. 

82.  Figure  25  is  the  slope  map  produced  wiien  the  smoothing  option 
was  selected  and  the  islands  were  removed  in  a subsequent  computer  run. 

Hunter-Liggett  Site  (1:25.000  Map  Scale) 

83.  The  steps  followed  to  produce  slope  maps  of  tiio  Hunter- 
Liggett  site  were  the  same  as  those  previously  used  to  produce  slope 
maps  of  the  Middle  East  site  (see  paragraph  71) • 

Site  elevation  data 

81*.  The  source  of  the  elevation  data  was  contour  lines  on 
1 :25 ,000-scale  U.  S.  Arm^  Topographic  Command  topographic  maps  identi- 
fied as  sheet  1755  IV  NE,  Series  V895S,  Edition  1-TPC,  Stony  Valley, 
California,  and  Sheet  1755  IflW,  Series  V895S,  Edition  1-TPC,  Jolon.  The 
contour  interval  used  on  these  maps  is  1*0  ft  with  supplementary  con- 
tours at  20  ft.  The  portion  of  the  maps  designated  in  this  study  as 
the  hunter- Liggett  site  is  shown  in  Figure  26.  Each  contour  line  with- 
in the  6.375-  by  2.825-km  site  was  digitized  and  processed  by  the  same 
technique  employed  for  the  Middle  East  site  to  generate  an  elevation 
grid  array.  Tlie  spacing  between  grid  points  was  set  this  time  at  a 
constant  25  m.  The  size  of  the  subsequent  array  was  256  by  111*  (NXRE 
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SLOPE  CLASS 
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0®-<5®  5®- <10®  10®- <30®  30®- <45®  45® -90® 

SLOPE  RANGE 

Figure  2U.  Slope  map  of  Middle  East  site  with  smoothing 
and  island  removal  options  not  requested  (map  produced 
with  the  WES  film  reader/writer 
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0®  -<5* 


5®- <10®  10®- <30®  30®- <45®  45® -90® 

SLOPE  RANGE 


Kif^ure  25.  Slope  map  of  Middle  East  site  witn  smoothing 
and  i.sland  removal  options  requested  (map  produced  with 
the  WES  film  reader/writer) 
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Figure  26.  Topographic  map  of  the  Hunter-Liggett  site  (map  scale  is  1:25,000) 


by  NYHE).  The  elevution  values  were  stored  as  integer  multiples  of 
decimetres,  and  tiie  DATASC  scale  constant  was  set  at  0.10. 

85.  Examination  of  the  perspective  plots  (Figure  27)  drawn  from 
the  Huiiter-Liggett  elevation  grid  array  reveals  the  same  terracing  ef- 
fect visible  in  the  Middle  East  perspective  plots  (Figure  19)  but  on  a 
much  reduced  scale.  This  is  true  because  the  Hunter-Liggett  elevation 
grid  array  was  constructed  using  a grid  spacing  constant  of  25  m,  which 
is  smaller  than  that  used  to  construct  the  Middle  East  elevation  grid 
array.  As  before,  however,  the  perspective  plots  show  that  the  Hunter- 
Liggett  elevation  grid  array  does  represent  to  a good  approximation  the 
three-dimensional  terrain  surface  of  the  Hunter-Liggett  site  and  is 
adequate  for  providing  input  topographic  data  to  SLOPEMAP. 

Selection  of  output  forms 

86.  The  class  ranges  listed  in  paragraph  7^  were  also  used  for 
this  site.  To  demonstrate  the  flexibility  of  SLOPfWAP,  the  following 
forms  of  output  were  chosen: 

a.  Swath  dump  of  classed  slope  values  only. 

b.  Magnetic  tape  of  classed  values  only. 

£.  No  punched  cards. 

Drum  plot . 

No  smoothing  options, 

£.  Slope  vector  components  on  magnetic  tape. 

The  control  values  necessary  to  run  each  of  these  selected  output 
options  were  again  obtained  from  Table  U.  For  these  selected  forms, 
the  output  option  card  (Figure  li4d)  was  set  up  as  follows: 

2, 2, 0,1, 0,1,  O.OOOOU 

87.  The  JCL,  data,  and  option  cards  that  were  necessary  to 
execute  the  example  with  a system  loadable  product  tape  are  shown 

in  Figure  28.  Note  that  for  illustrating  the  calculation  units,  the 
previously  selected  class  ranges  in  degrees  are  now  expressed  as 
tangents  of  the  angles.  The  slope  class  boundaries  and  areas  would 
have  been  the  same  if  degrees  had  been  used  instead  of  tangent  values. 
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Figure  28.  Job  control  language,  data,  and  option  cards  used 
to  produce  the  slope  map  shown  in  Figure  29  of  the  Hunter- 

Liggett  site 
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Slope  map  pictures 

8fl.  The  cla3sed  value  tape  was  used  again  as  an  input  tape  to 
the  WES  film  reader/writer  to  produce  a slope  map  picture  (see  para- 
grapli  8l ) . The  slope  map  produced  from  this  process  of  the  Hunter- 
Liggett  site  is  shown  in  Figure  29.  Isolated  classed  slope  value  islands 
and  terracing  effects  that  tend  to  complicate  the  appearance  of  the  slope 
map  are  less  evident  than  they  were  in  the  Middle  East  site  (Figure  2k). 
However,  after  SLOPEMAl  was  rerun  with  both  smoothing  options,  the  gen- 
erated slope  map  (Figure  30)  became  less  complicated. 

Comparison  of  Slope  Maps  Constructed  by 
Manual  and  Computer  Methods 


89.  To  compare  slope  maps  constructed  by  a manual  method  with 
those  constructed  by  the  computer  method  developed  in  this  study,  both 
manually  and  computer-produced  slope  maps  of  the  Middle  East  and  Hunter- 
Liggett  sites  were  constructed  at  the  same  map  scale. 

Manual  procedures 

90.  Manually  prepared  slope  maps  of  both  the  Middle  East  and 
Hunter-Liggett  sites  were  constructed  using  the  procedure  discussed  in 
Reference  11.  Construction  was  begun  by  selecting  the  map  scale  and 
contour  interval  on  templates  (Figure  31)  that  coincided  with  that  of 
the  maps.  The  templates  were  transparent  so  that  contour  lines  could  be 
readily  seen  through  them.  This  was  necessary  since  the  technique  of 
ascertaining  slope  classes  at  everv  location  required  overlaying  the 
template  on  the  map  and  moving  it  freely  over  the  contour  lines.  The 
slope  classes  were  determined  by  finding  the  set  of  slope  class  circles 
or  lines  (sequence  of  vertical  lines  immediately  to  the  left  of  the 
slope  class  circles)  that  bracketed  the  separation  between  two  contour 
lines.  For  example,  if  at  a given  location  on  the  topographic  sheet  the 
distance  between  a pair  of  contour  lines  was  greater  than  the  inner 
circle  of  a given  slope  class  range  but  less  than  the  outer  circle  for 
that  range,  the  map  location  was  classed  as  that  slope  class  range  and 
marked  accordingly.  This  process  was  repeated  many  times  over  the 
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I.  Slope  map  of  the  Hunter-Liggett  site  with  smoothing  and  island  removal 
options  requested  (map  produced  with  the  WES  film  reader/writer ) 
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Figure  31.  Template  used  to  manually  construct  slope  maps 
from  the  topographic  data  in  Figures  l8  and  26 

entire  topographic  map  until  the  slope  conditions  had  all  been  determined. 
However,  this  process  was  very  subjective  at  times,  and  some  error  in 
interpretation  was  expected.  The  procedure  using  vertical  lines  is,  in 
principle,  the  same  as  for  circles.  Note  that  the  reference  line  of 
each  set  is  always  the  leftmost  line  of  the  set. 

Comparison  of  manually  and 
computer-produced  slope  maps 

91.  The  computer-generated  areas  corresponding  to  slope  classes 
A , B , and  C for  both  tho'  Middle  East  and  Hunter-Liggett  sites 
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are  presented  alonj!;  witti  their  manually  produced  counterparts  in  Vip^- 
ures  i?~31  • Classes  D 'ind  K were  not  compared  in  this  study  since 
as  a sum,  they  j-epresented  less  than  0.1  percent  of  the  total  areas 
mapped,  and  as  individual  contiguous  patches,  they  represented  areas  no 
larger  than  0.01  percent  of  tiie  total  areas.  Valid  comparison,  it  was 
felt,  co’.ila  not  be  made  for  such  small  areas,  since  mapping  oversights 
due  to  tne  subjective  pr..'cessos  inherent  in  the  manual  methods  of  con- 
structing slope  maps  caust  i frequent  omissions  of  small  areas.  In  view 
of  these  Svibjective  processes,  the  spatial  placements  of  the  different 
classed  areas  by  the  two  methods  were  in  remarkably  close  agreement.  The 
major  disagreement  in  these  comparisoiis  originates  from  the  artificial 
terracing  in  the  elevation  grid  array.  This  effect  is  most  readily 
observed  in  comparison  of  slope  classes  A and  B . 

92.  Comparison  of  slope  class  areas  resulting  from  the  manual 
and  computer  construction  methods  indicates  that  the  computer  methods 
generated  comparable  results.  The  areas  of  classes  A , B , and  C 
as  determined  by  each  of  the  two  construction  methods  are  presented  in 
Table  6.  Units  of  area  are  expressed  as  a percentage  ol'  the  total 
mapped  area,  and  the  criterion  for  comparison  is  percent  deviation 
between  the  two  methods  of  construction.  Percent  deviation  is  defined 
as  |A  -A1/(A  +A)x  100%  , where  A and  A are  the  areas  of 
the  computer  and  manually  produced  slope  classes,  respectively.  Table  6 
shows  that  the  areal  calculations  agree  very  well  except  within  class  B. 
This  is  not  surprising,  however,  since  it  was  within  this  class  that  the 
terracing  effect  was  most  prevalent. 

Computer  Core  Space.  Time,  and  Cost  Requirements 
for  Executing  SLOPEMAi’ 

93.  Estimates  of  computer  core  space,  running  times,  and  costs 
required  for  constructing  slope  maps  using  the  automated  procedure  are 
presented  below.  The  costs  of  constructing  slope  maps  by  manual  and 
computer  methods  are  also  compared. 
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a . Manual 


b.  Computer 


Figure  33.  Comparison  of  slope  class  B (5-‘^10  deg,  in  black) 
in  the  Middle  East  site  mapped  by  manual  and  computer  methods 
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a . Manual 


b.  Computer 

Fig\ire  3^.  Comparison  of  slope  class  C (l0-<30  deg,  in  black) 
in  the  Middle  East  site  mapped  by  manual  and  computer  methods 


a . Manual 


b.  Computer 


Figure  35*  Comparison  of  slope  class  A (0-<5  deg, 
in  black)  in  the  Hunter- Liggett  site  mapped  by  manual 
and  computer  methods 


b.  Computer 


Figure  36.  Comparison  of  slope  class  B (5-<10  deg, 
in  black)  in  the  Hunter-Liggett  site  mapped  by  manual 
and  computer  methods 
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b.  Computer 


Figure  37.  Comparison  of  slope  class  C (l0-<30  deg, 
in  black)  in  the  Hunter-Liggett  site  mapped  by  manual 
and  computer  methods 
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Core  space  requirements 

9^*.  The  necessary  core  space  requirements  should  be  known  prior 
to  running  SLOPKMAI'  so  that  the  correct  specifications  can  be  made  on 
the  $:LIMITS:  and  $:FILE;  JCL  cards  (Figure  l6).  Compilation  and 
load  requirements  are  specified  on  the  first  and  second  $:LIMITS: 
cards  (lines  50  and  110,  Figure  l6),  respectively.  The  requirements  for 
the  main  program  and  each  subroutine  can  be  determined  from  Table  7 for 
the  WES  Honeywell  G-635  computer.  The  necessary  core  space  required  to 
compile  is  determined  by  the  largest  of  the  routines  of  the  program. 

Since  the  main  program  is  always  required  and  is  the  largest  routine, 
the  limits  to  compile  should  be  set  large  enough  to  accommodate  37K 
words  of  space  (a  word  on  the  Honeywell  G-635  system  consists  of  a bit 
string  of  36  binary  bits).  Limits  for  program  loading  should  be  set  to 
accommodate  the  accumulative  load  requirements  of  all  the  subroutines. 

For  example,  a limit  equal  to  or  greater  than  3^K  words  should  be 
specified  for  loading  all  of  program  SLOPEMAP. 

95.  The  random  access  file  space  requirements  for  file  codes  30 
and  31  (lines  250  and  260,  Figure  l6)  can  be  determined  from  the  total 
number  of  grid  locations  in  the  input  elevation  grid  array.  For  an 
array  with  NXRE  columns  and  NYRE  rows,  the  total  requirement  is  the 
product  NXRE  x NYRE  expressed  in  terms  of  links  (a  link  is  38^40  words). 
For  file  32  (line  270,  Figure  l6),  the  file  space  required  is  twice  the 
product  of  NXRE  and  NYRE  expressed  in  terms  of  links. 

Processing  times 

96.  The  time  required  to  compile  SLOPEMAP  is  nominally  35  sec  on 
the  WES  Honeywell  G-635  computer.  This  time  can  double  if  reference 
maps  and  stoi-age  maps  are  requested  as  additional  output  forms  to  the 
compilation  reports.  This  fact  is  of  value  only  if  the  program  is 
modi fied. 

97.  Processing  times,  once  the  program  begins  slope  calculations, 
depend  on  the  number  and  type  of  output  options  selected,  as  well  as  on 
the  total  number  of  grid  locations  in  the  input  grid  array.  Runs  were 
made  with  five  grid  arrays  whose  dimensions  varied  from  100  by  100  up 

to  500  by  500  to  establish  an  approximate  dependence  of  processing  times 
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on  size  of  the  grid  array.  The  runs  were  made  with  options  to  output  a 
classed  value  swath  dump  and  a magnetic  tape  of  classed  values.  Two 
runs  were  made  for  each  of  the  five  grid  arrays,  one  with  both  smooth- 
ing options  requested  and  one  without.  The  results  of  these  runs  are 
plotted  in  Figure  38  and  should  provide  a basis  for  estimating  process- 
ing times. 

Processing  costs 

98.  The  cost  estimates  for  running  SLOPEMAP  are,  of  course, 
computer  and  installation  dependent.  Since  all  of  the  development  work 
was  accomplished  with  the  WES  Honeywell  computer  system,  only  cost 
estimates  based  on  the  1976  WES-system  charges  are  presented.  Compila- 
tion and  loading  charges  for  SLOPEMAP  totaled  approximately  $U.60  with- 
out requesting  reference  and  storage  maps  and  $9.00  if  maps  were  re- 
quested. These  charges  were  not  incurred  if  the  program  was  executed 
each  time  from  a system  loadable  product  tape  (see  paragraph  68). 

99.  Processing  costs  as  a function  of  input  grid  array  size  were 
derived  from  the  same  set  of  computer  runs  used  to  develop  Figure  38; 
that  is,  the  costs  as  a function  of  input  grid  size  were  determined  for 
outputting  a classed  value  swath  dump  and  a magnetic  tape  of  classed 
values  for  both  smoothing  and  no  smoothing  requests.  These  costs  are 
presented  in  Figure  39  as  a function  of  the  total  number  of  elevation 


grid  locations.  Computer  processing  costs  per  unit  area  were  also  de- 
termined from  these  data  and  are  presented  in  Figure  hO  as  a function 
of  the  grid  spacing  constants.  Cost  calculations  using  Figure  hO  must 
include  $2.00  for  peripheral  computer  hardware  (tapes  and  high-speed 
printer).  No  compilation  charges  are  included  in  these  curves  since 
the  rxins  were  all  made  with  a system  loadable  product  tape. 

100.  llie  additional  costs  of  running  SLOPEMAP  with  other  output 
options  were  determined  from  the  example  runs  of  the  Middle  East  site 
that  produced  the  slope  map  shown  in  Figure  2^4 . The  total  cost  of 
running  SLOPEMAP  for  this  Job  was  about  $26.00.  The  greater  part  of  the 
increase  in  cost  over  an  estimate  obtained  from  Figure  39  was  from  the 
punching  of  cards  that  cost  about  $7.00.  The  remaining  cost  was  $5.00 
for  the  CRT  plot,  $14.00  for  the  slope  vector  tape,  and  $2.00  for  the 
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PROCESSING  TIME.  SEC 


TOTAL  NUMBER  OF  ELEVATION  GRID  LOCATIONS 


Figure  38.  Functional  dependence  of  SLOPEMAP  computer  processing 
time  (with  and  without  smoothing  options)  on  the  input  elevation 

grid  array  size 
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Figure  39.  SLOPEMAP  computer  processing  costs  (with  and  without 
smoothing  options)  as  a fxinction  of  the  total  number  of  input 

elevation  grid  locations 
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COMPUTER  PROCESSING  COST.  DO 


Figure  Uo.  SLOPEMAP  computer  processing  costs  (with  and  without  smooth- 
ing options)  as  a function  of  the  grid  spacing  constant 
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calculated  slope  value  swath  dump.  Running  this  same  job  with  both 
smoothing  options  increased  the  total  cost  about  $U.OO. 

Comparison  of  costs 

101.  The  costs  for  constructing  slope  maps  for  this  study  auto- 
matically and  manually  can  be  compared  by  using  the  two  study  sites. 
Middle  East  and  Hunter-Liggett  (see  paragraphs  71-88).  Since  the  two 
sites  were  approximately  the  same  size  in  terms  of  total  grid  points, 
the  costs  for  automatically  constructing  slope  maps  for  them  were  about 
equal.  Therefore,  the  following  costs  presented  are  representative 

of  both  sites. 

102.  The  total  cost  of  the  manually  constructed  slope  maps  with 
five  delineated  classes  was  approximately  $150.00  (based  on  one  man-day 
of  labor).  The  cost  of  changing  slope  class  ranges  and  constructing 
yet  another  slope  map  with  five  classes  was  also  $150.00,  even  for  the 
same  topographic  data.  This  is  an  important  point  to  make  for  compari- 
son, since  repeating  computer  procedures  to  perform  reclassification 
represents  little  additional  cost.  This  is  the  real  payoff,  since  the 
selection  of  class  ranges  is  very  subjective  and  the  slope  categories 
would  more  than  likely  be  different  for  different  types  of  project 
assignments  that  required  such  data. 

103.  The  total  cost  of  the  slope  maps  constructed  by  automated 
means  in  this  study  was  the  siim  of  the  costs  of  four  production  steps, 
as  follows; 

Production  Step 

1.  Digitization  of  contour  lines  of 
input  topographic  map 

2.  Execution  of  computer  program  to 
generate  elevation  grid  array 

3.  Execution  of  SLOPEMAP  to  produce 
classed  value  swath  dump  and 
magnetic  tape 

h.  Production  of  a black-and-white 
transparency  on  reader/writer 

The  sum  of  these  production  steps,  which  represent  the  cost  of  producing 
the  first  computer  slope  map,  is  $120.00.  Note  that  succeeding  slope 
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Approximate  Cost 
$75.00 

25.00 

12.00 

8.00 


maps  produced  from  the  same  input  elevation  grid  array  require  only  an 
additional  $20.00  each.  Hence,  reclassification  of  given  areas  to  ob- 
tain slope  data  for  new  baseline  data  requirements  by  computer  methods 
becomes  a very  inexpensive  process  when  compared  with  the  costs  of 
manually  performing  the  same  task.  Further  savings  can  be  realized 
during  subsequent  reclassification  runs,  if  the  smoothing  options  are 
needed,  by  selecting  as  an  output  product  a vector  component  tape  during 
the  initial  run  of  SLOPfMAF  and  then  using  this  tape  as  the  input  topo- 
graphic data  source  during  the  subsequent  runs.  This  avoids  the  neces- 
sity of  having  to  rerun  with  the  smoothing  options  during  each  reclas- 
sification and  can  represent  a considerable  savings  (Figure  39)  For 
large  arrays. 


PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 
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♦ 


Conclusions 


lOU.  Conclusions  drawn  from  this  study  are: 

a.  The  automated  procedure  (SLOPEMAP)  developed  in  this 
study  provides  the  Army  and  Research  Community  with  a 
cost-effective  way  of  generating  reliable  and  accurate 
baseline  data  on  slope.  The  scope  of  quantitative 
methods  for  assessing  the  environmental  impact  of  mili- 
tary installation  operations  that  require  slope  as  a 
baseline  data  element  have  increased  steadily.  Digital 
data  and  maps  of  up  to  36  user-selected  class  values  of 
slope  can  be  obtained  with  SLOPEMAP  making  it  very  flex- 
ible for  satisfying  specific  data  requirements. 

b.  The  initial  cost  of  constructing  slope  maps  using  SLOPE- 
MAP  is  comparable  to  manual  methods;  however,  the  cost  of 
constructing  reclassified  slope  maps  to  support  new  base- 
line data  requirements  using  previously  generated  topo- 
graphic input  data  can  be  from  five  to  six  times  less 
using  SLOPEMAP  than  the  cost  of  manual  methods  (see 
paragraphs  102  and  103). 

c_.  A vector  component  tape  should  be  output  during  the 
initial  run  of  SLOPEMAP  if  smoothing  options  are 
selected  and  if  subsequent  reclassification  runs  are  to 
be  made.  Reclassification  runs  should  use  this  tape  as 
the  input  topographic  data  source  in  order  to  reduce  the 
time  and  cost  of  execution.  This  source  of  input 
prevents  the  computer  from  having  to  duplicate  numerous 
calculations  (see  paragraph  38). 

d.  CRT  plots  should  be  used  during  the  early  stages  of 
slope -map  construction.  After  output  needs  have  been 
established,  the  more  accurate  and  more  costly  drum  plot 
option  can  be  used  (see  paragraph  57). 

Recommendations 

105.  Recommendations  that  would  increase  the  accuracy  and  utility 
of  SLOPEMAP  are: 

a.  More  accurate,  reliable,  and  cost-effective  techniques 
should  be  developed  for  producing  elevation  grid  arrays 
from  primary  topographic  data  sources. 

SLOPEMAP  should  be  made  directly  compatible  with  the 
DMATC  topographic  data  tapes  (see  paragraph  k2). 
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Table  1 


farllal  Herlvatlves  of  f(x.y)  Correapondlnn  to  the 
Algorlttis  for  Calculatltm  Slope 


Vector 

Plane 

I’urface 

^r,Ox 

»f/Jy 

1 

Vector  Algorithm* 

(z^  - Zq)/D 

0 

2 

— 

— 

^*2  ■ *0^^^ 

(Zg  . Zo)/2D 

3 

— 

— 

0 

(Z3  - Zq)/D 

U 

— 

— 

(z,^  - Zp)/20 

5 

— 

— 

(Zq  - Zj)/D 

0 

6 

— 

— 

(Zq  - Zf)l2Q 

(Zq  - Zg)/2D 

7 

— 

— 

0 

(Zq  - z^)/D 

8 

— 

— 

(Zg  - Zq)/2D 

(Zq  - Zg)/2D 

1 

• • 

Plane  Algorithm 

(z,  - Zq)/D 

(Z3  - Zgj/D 

— 

2 

— 

(Zq  - *5)/^ 

(Z3  - Zg)/D 

— 

3 

— 

(Zq  - *5)/^ 

(Zg  - z^)/D 

— 

— 

(z^  - 

(Zq  - z^)/D 

S 

Surface  Algorithm 

(z^  - z^)/2D 

(z^  - z^)/2D 

S' 

[ (*2  - Zg)  - (Zj^  - Zg)]/1*D 

((Z2  - Zg)  ♦ (Zjj  - Zq))/**D 

' 

* D Is  the  grid  spacing  constant,  and  z represents  elevation  values  at  the 
different  grid  locations  as  designed  for  the  vector,  plane,  and  surface 
algorithms  in  Figures  3,  **,  and  5,  respectively. 


Table  2 


Comparisons*  of  the  Vector.  Plane,  and  Surface  /Ugorlthm 
Slope  Values  with  the  IJcact  Values 


Percent  Probability  that  the  Difference  Between 
the  Algorithmic  and  Exact  Slopes  Fell  Within  a 
Specific  Slope  Difference  Range 

Slope  Difference 
Hange.  Deg 

Vector 

Algor  ititm 

Plane 

Algorithm 

Surface 

Algorithm 

0 to  <1 

25.3 

21.9 

93.8 

1 to  <2 

1*3.2 

1*0.1* 

6.2 

2 tc  <U 

19.1 

25.3 

0.0 

U to  <6 

10.9 

10.9 

0.0 

6 to  <11 

1.5 

1.5 

0.0 

• All  comparisons  were  made  with  eurrays  containing  100  by  100  grid 
points. 
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Table  3 

Index  of  DMATC  Large-Scale  Digital  Terrain  Elevation  Data 
on  Military  Installations  in  the  United  States 


o* 


Area 

Scale 

Coordinates 

AZ, 

Fort  Huachuca 

1:25,000 

Lat  N 

31®22'30” 

31°l45’00" 

Long  W 

110°07'30" 

- 

110° 37 ’30" 

CA, 

Camp  Pendleton 

1:214,000 

Lat  N 

33®12'00" 

33° 36 *00" 

Long  W 

117“12'00" 

- 

117°36'00" 

CA, 

Fort  Irwin 

1:U9,212 

Lat  N 

35007’ 30” 

35° 37 ’30" 

Long  W ll6®17'30" 

- 

ll6°57'30" 

CA, 

Fullerton 

1:149,212 

Lat  N 

33052*  30" 

314°  26 ’00" 

Long  W 

117°25'00" 

- 

117° 57 ’00" 

CA, 

Fort  Hunter- 

1:50,000 

Lat  N 

35045 '00" 

36°00’00" 

Liggett 

Long  W 

121° 00 ’00" 

- 

121°15’00" 

o 

> 

Fort  Hunter- 

1:25,000 

Lat  N 

35® 52 ’30" 

36°00'00" 

Liggett 

Long  W 

121° 07 ’30" 

- 

121°15'00" 

KA, 

Fort  Riley 

1:50,000 

Lat  N 

39°10'00" 

39°20' 00" 

Long  W 

96°l45'00" 

- 

97°00’00" 

NC, 

Camp  Leje’one 

1:214,000 

Lat  N 

314°  30 ’00" 

3l4°l48’00" 

Long  W 

77°06’00" 

- 

77° 30 ’00" 

NC, 

Camp  Lejexuie 

1:50,000 

Lat  N 

314°  30  *00" 

3l4°l48’00" 

Long  W 

77° 00 ’00" 

- 

77° 30 ’00" 

OK, 

Fort  Sill 

1:50,000 

Lat  N 

3l4°32’00" 

3l4°l49'56" 

Long  W 

98°15'00" 

- 

98°l47'15" 

TX, 

Fort  Hood 

1:50,000 

Lat  N 

31°00'00" 

31°30'00" 

Long  W 

97° 30 ’00" 

- 

98°00'00" 

< 

> 

Quant Ico 

1:214,000 

Lat  N 

38"2l4'00" 

* 

38°l42'00" 

Long  W 

77°12’00" 

- 

77° 36 *00" 

Table  1* 


Control  Values  for  Output  Option  Card 


Card 

Control 

Output  Option 

Field* 

Value 

Output  Description 

Swath  dump 

1 

0 

Ho  output 

1 

Calculated  values  only 

2 

Classed  values  only 

3 

Calculated  and  classed  values 

Magnetic  tape 

2 

0 

No  output 

1 

Calculated  values  only 

2 

Classed  values  only 

3 

Calculated  and  classed  values 

Punched  cards 

3 

0 

No  output 

1 

Calculated  values  only 

2 

Classed  values  only 

3 

Calculated  and  classed  values 

Drum  or  CRT  plot 

0 

No  output 

1 

Drum  plot  only 

2 

CRT  plot  only 

Smoothing  options 

5 

0 

No  processing 

1 

Smoothing  only 

2 

Island  removal  only 

3 

Smoothing  followed  by  island 

removal 

Vector  components 

6 

0 

No  output 

1 

Output  slope  vector  components 

onto  magnetic  tape 

Drum  plot  scale 

7 

Drum  plot  scale  value  in  deci- 

mal form  (e.g. 

1:25000  -«•  O.OOOOU) 

i 


* Card  fields  are  defined  by  the  field  delimiters  (e.g.  blanks  and 
commas).  The  length  of  each  field  depends  solely  on  the  placement  of 
the  delimiter  on  the  card. 

I 

» 


Table  5 


File  Codes  and 

Descriptions  of  the 

Input  and  Output 

Files  Used  in  SLOi'iWAi’ 

Function 

File  Code 

Tyi.e  File* 

Description  of  Contents 

Input 

10 

Sequential 

Elevation  data 

Input 

11 

Sequential 

Slope  vector  data 

Output 

20 

Sequential 

Calculated  slope  values 

Output 

21 

Sequential 

Classed  slope  values 

Output 

22 

Sequential 

Slo|>e  vector  components 

Output 

23 

Sequential 

Drum  plot 

Output 

2t* 

Sequential 

CRT  plot 

Scratch 

30 

Random 

Calculated  slope  values 

Scratch 

31 

Random 

Classed  slope  values 

Scratch 

32 

Random 

Slope  vector  components 

» The  sequential  files  most  frequently  used  during  SLOPEMAP  development 
were  magnetic  tape  and  pennanent  disc  files  (PRMFL). 


Table  6 


Comparison  of  Computer  and  Manually  Classed 
Sloi>e  Class  Area 


Area**  of  Computer  Area**  of  Manually  Percent 

Class*  Classed  Slope  Class  Classed  Clope  Class  rievlatlon 

Middle  tost  Site 


A 

51.0 

55.1 

%9 

B 

29.1 

23.2 

11.3 

C 

19.8 

21.1* 

3.9 

Hunter-Llaaett  Site 

A 

57.1 

65.1* 

00 

NO 

B 

20.9 

11*. 8 

17.1 

C 

22.0 

19.8 

5.3 

• Classes  D and  E represented  less  than  0.1  percent  of  the  total 
areas  classed  and  thus  were  not  compared. 

**  Areas  are  expressed  in  percentage  of  total  area  mapped. 


Table  7 


SLOPKMAP  Computer  Core  Space  Requlrementa 


Subroutine 

Name 

Subroutine  Function 

Compilation 

Requirements* 

Load 

Requirements' 

MAIN 

Calculates  slopes 

37K 

18K 

SMOOTH 

Averages  slope  vector  components 

28k 

2K 

ISLAND 

Removes  classed  value  islands 

27K 

IK 

CALi'WH 

Prints  calculated  slope  values 

27K 

2K 

CLSSWH 

Prints  classed  slope  values 

32K 

2K 

CRTPLT 

Produces  a CRT  plot 

29K 

UK 

DRMPLT 

Produces  a DRUM  plot 

28k 

2K 

CRDCAL 

Punches  a card  deck  of 
calculated  slope  values 

26k 

1/2K 

CRDCLS 

Punches  a card  deck  of 
classed  slope  values 

26k 

1/2K 

TAPCAL 

Produces  a magnetic  ta'pe  of 
calculated  values 

25K 

1/2K 

TAPCLS 

Produces  a magnetic  tape  of 
classed  values 

25K 

1/2K 

TAPVEC 

Produces  a magnetic  tape  of 
slope  vector  components 

25K 

IK 

• Units  for  compilation  and  load  requirements  are  given  in  terms  of 
Honeywell  0-635  system  words.  Each  word  contains  a string  of  36 
binary  bits. 
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In  accordance  with  letter  from  OAEN-ROC,  DAEN-ASl  dated 
22  July  1977,  Subject:  Facsinlle  Catalog  Cards  for 
Laboratory  Technical  Publications,  a facsinlle  catalog 
card  In  Library  of  Congress  NARC  fomat  is  reproduced 
below. 


Struve,  Horton 

An  autouted  procedure  for  elope  nap  construction  / by 
Horton  Struve.  Vicksburg,  Mias.  : U.  S.  Uaterwaya  Eaperlnent 
Station.  1977. 

2 V.  : 111.  ; 27  cn.  (Technical  report  - U.  S.  Amy  Engineer 
Waterways  Eaperlnent  Station  ; M-77-3) 

Prepared  for  Office,  Chief  of  Engineers,  U.  S.  Amy,  Wash- 
ington, D.  C. , under  Project  4A1S212IA896,  Task  01,  Work  Unit 
006.  Project  1E86S803H710. 

Includes  bibliography. 

Contents:  v.l.  Description  and  instructions  for  use  of  the 
automated  procedure.  (Distribution  unllnited)  .— v.  2.  Listing 
and  glossary  for  program  SLOPEMAP.  (Distribution  limited  to 
U.  S.  Government  agencies  only). 


